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FOREWORD 


This  report  deals  with  the  prediction  of  explosion  injury  to  fish  with 
swim  bladders  and  is  part  of  a  continuing  study  of  the  effects  of 
underwater  explosions  on  marine  life.  Swim  bladder  fish  are  parti¬ 
cularly  vulnerable  to  explosions,  and  this  group  includes  the 
majority  of  fish  with  sports  and  commercial  value.  This  study  will 
result  in  an  improved  capability  to  predict  such  effects,  and  will  be 
useful  in  connection  with  the  testing  of  new  explosives  and  of  warheads 
at  sea. 

This  study  is  part  of  the  ordnance  pollution  abatement  program  of  the 
Naval  Sea  Systems  Command  and  was  supported  by  SEA  TASK  SSL  55001/ 
19373. 

The  author  is  indebted  to  Ermine  A.  Christian  and  George  A.  Young  for 
many  valuable  suggestions  during  the  course  of  this  work. 
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bladder  to  the  incident  explosion  pressure  wave.  How  these  injuries  come  about- 
whether  they  are  due  to  displacements  associated  with  the  motion  or  to  the  radiated 
pressure  field-is  not  known. 
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DYNAMICAL  MODEL  FOR  EXPLOSION  INJURY  TO  FISH 

1  INTRODUCTION  AND  SUMMARY 

This  report  describes  a  method  for  calculating  the  kill  probability 
for  bladder  fish  subjected  to  an  underwater  explosion.  The  method 
consists  of  an  approximate  calculation  for  the  extreme  values  of  com¬ 
pression  and  extension  of  the  fishes'  gas-filled  swim  bladder  in 
response  to  the  explosion  pressure  wave.  The  calculations  are  made 
for  the  damped  radial  oscillations  of  a  spherical  air  bubble  in  water. 
The  kill  probability  is  then  calculated  as  an  experimentally  determined 
function  of  the  ratio  of  maximum  to  minimum  radius  during  the  oscilla¬ 
tory  response.  For  each  species  of  fish  the  effective  swim  bladder 
radius  needed  for  the  calculations  is  determined  by  trial  and  error 
thru  correlating  the  observed  injuries  with  the  calculated  ratio  of 
maximum  to  minimum  radius. 

The  method  is  used  to  describe  dissection  results  from  1500  Spot 
and  White  Perch  caged  at  depths  from  5  to  100  feet  and  subjected  to 
pentolite  explosions  of  from  1  to  70  lbs  submerged  at  depths  ranging 
from  from  5  to  70  feet.1,? 

These  results  are  then  used  to  calculate  several  examples  of 
kill-probability  and  fish-kill  distributions  for  a  particular 
explosion  geometry  and  several  species  and  sizes  of  fish.  Details 
of  the  computations  are  given  in  the  appendices.  Appendix  C  gives 
an  approximate  method  for  calculating  the  pressure  signature  from  an 
underwater  explosion  subject  to  surface  effects  including  cavitation. 


1.  Gaspin,  Joel  B.,  "Experimental  Investigations  of  the  Effects  of  Underwater 
Explosions  on  Swim  Bladder  Fish,  I:  1973  Chesapeake  Bay  Tests",  NSWC/WOl.  TR  75-58, 
1975. 

2.  Gaspin,  J.  B.,  Wiley,  M.  L.,  and  Peters,  G.  B.,  "Experimental  Investigations  of 
the  Effects  of  Underwater  Explosions  on  Swim  Bladder  Fish,  II:  1975  Chesapeake  Bay 
Tests",  NSWC/WOL  TR  76-61,  1976. 
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2 .  DYNAMICAL  MODEL 

2.1  THE  SWIM  BLADDER.  Figure  2.1.1  shows  a  life-size  tracing  of 
a  Spot  (Leiostomus  xanthurus  Lacep'ede)  which  depicts  the  swim  bladder 
and  the  surrounding  internal  organs.  Figure  2.1.2  shows  a  swim  bladder 
removed  from  a  Spot.  In  our  work  we  depicted  the  motion  of  this  organ 
by  that  calculated  for  a  spherical  bubble  of  air.  The  size  of  this 
"equivalent  spherical  bubble"  was  first  estimated  by  equating  its 
volume  to  that  of  the  fishes'  swim  bladder  —  then  later  we  adjusted 
the  size  of  the  equivalent  bubble  to  optimize  the  correlation  between 
the  calculated  motion  and  the  observed  injuries. 

Figure  2.1.3  shows  a  life-size  tracing  of  a  White  Perch 
(Mo rone  americana  Gmelin) ,  the  other  species  of  fish  used  in  the  work 
reported  here.  The  inset  shows  the  shape  of  the  forward  portion  of 
the  swim  bladder  in-place  in  the  fish. 

2.2  PRELIMINARY  CALCULATIONS  OF  BLADDER  OSCILLATION.  We  began 
this  study  by  making  order  of  magnitude  calculations  which  approxi¬ 
mated  the  explosion  pressure  field  with  a  step  increase  followed  by  a 
step  decrease  in  pressure.  The  fishes'  swim  bladder  was  approximated 
by  a  spherical  bubble  of  air  in  an  infinite  body  of  water  —  and,  we 
calculated  its  undamped  oscillatory  response.  These  order  of  magnitude 
calculations  showed  that  the  oscillatory  response  of  the  swim  bladder 
was  a  likely  source  of  the  fishes'  injuries.  They  also  pointed  out  a 
strong  resonance  which  occurs  when  surface  cut-off  (the  arrival  of  the 
rarefaction  wave  reflected  from  the  water  surface)  happens  at  the 
instant  of  maximum  bladder  compression.  This  work  was  described  in  an 
internal  Center  report  which  is  reproduced  as  Appendix  A  of  this  report. 

2.3  CURRENT  METHOD  FOR  CALCULATING  BLADDER  OSCILLATION.  The  method 
developed  for  calculation  of  the  oscillatory  response  to  the  changing 
pressure  field  generated  by  an  underwater  explosion  is  described  in 
detail  in  Appendix  B.  Here,  we  attempt  to  give  the  reader  a  feel  for 
this  computation  and  its  approximate  nature.  Figure  2.3.1  shows  the 
idealized  pressure  signature  measured  on  an  underwater  explosion  test 
together  with  the  calculated  pressure  inside  a  fish's  swim  bladder 


This  is  a  life-size  body  tracing  of  a  frozen  fish  which  has  been  sawed  in  half. 


PRESSURE  IN  BLADDER 


BLADDER  PRESSURE  AND  SIZE  AS  A  FUNCTION  OF  TIME 
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as  it  responds  to  this  applied  pressure.  Below  the  pressure  traces  are 
corresponding  sizes  calculated  for  the  spherical  air  bubble  used  to 
represent  the  fish's  swim  bladder. 


This  fluid  motion  can  be  likened  to  a  mass  attached  between 
two  springs,  the  end  of  one  being  fixed,  the  end  of  the  other  movable 


MASS 


The  first  spring  is  analogous  to  the  internal  qas  pressure;  the 
second  to  the  changing  external  water  pressure.  Let  us  imagine  the 
mass  resting  on  a  frictionless  block  of  ice.  It  is  at  rest  in  a 
position  of  equilibrium.  If  at  t=0  we  quickly  displace  the  movable 
spring  to  a  new  position,  the  mass  will  oscillate  about  a  new  position 
of  equilibrium  determined  by  the  displacement  of  the  end  of  the  spring 
The  mass  will  oscillate  about  this  new  equilibrium  position  with  an 
amplitude  equal  to  the  displacement  of  this  equilibrium  position  from 
the  initial  at-rest  position  of  the  mass. 


This  initial  at-rest  position  now  becomes  an  extremum  of  the 
oscillatory  motion.  Since,  at  each  extremum  of  the  motion  the  mass 
is  instantaneously  at  rest,  a  subsequent  jump  displacement  of  the  end 
of  the  spring  occurring  at  an  extremum  simply  shifts  the  motion  to  an 
oscillation  about  a  new  equilibrium  position  -  and  the  new  amplitude 
is  just  the  distance  of  the  mass  at  the  time  of  the  jump  from  this 
new  equilibrium  position. 


Likewise,  for  the  bubble,  step  changes  in  the  outside  water 
pressure  occurring  at  half-period  intervals  simply  change  the  equili 
brium  pressure  (equal  to  the  outside  pressure)  of  the  oscillating 
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bubble  flow  -  and  the  new  amplitude  is  just  the  instantaneous 
excursion  of  the  bubble  pressure  from  the  changed  outside  pressure. 

We  made  use  of  this  property  of  the  motion  to  calculate  the  motion  as 
the  response  to  a  series  of  pressure  steps  as  shown  in  Figure  2.3.2. 
(Appendix  B  describes  the  way  we  estimated  successive  half-periods 
in  order  to  construct  the  approximating  step  wave.) 

The  transition  from  the  positive  to  the  negative  pressure 
phase,  i.e.,  surface  cut-off,  must  be  treated  as  a  special  case,  since 
in  general,  this  will  not  occur  at  an  extremum  of  the  motion.  If 
cut-off  occurs  at  a  compression,  this  compression  is  then  the  compres¬ 
sion  of  the  next  half  cycle  also  -  and  the  subsequent  expansion  is 
the  greatest  possible.  If  cut-off  occurs  at  an  expansion  and  the  new 
outside  pressure  is  lower  than  that  inside,  this  expansion  becomes 
the  compression  of  the  next  half  cycle,  there  occurs  a  180  degree 
phase  shift,  and  the  subsequent  expansion  is  the  smallest  possible. 

In  general,  cut-off  occurs  somewhere  in  between  and  there  is  an 
intermediate  phase  shift  and  subsequent  expansion  or  compression. 

In  such  cases  we  calculated  this  next  half  cycle  at  lowered  pressure 
as  starting  at  the  time  of  the  previous  extremum  but  at  a  ficticious 
amplitude  interpolated  as  a  function  of  the  time  at  which  cut-off 
occurred.  (The  details  of  this  procedure  are  described  in  Appendix  B.) 

Finally,  at  the  end  of  the  negative  phase,  the  return  to 
ambient  pressure  can  cut  short  or  even  prevent  the  final  expansion. 

This  was  handled  in  an  analogous  manner  to  surface  cut-off  and  is 
also  described  in  Appendix  B. 

To  approximate  the  dissipation  of  energy  which  in  nature  must 
occur  during  the  oscillation,  we  arbitrarily  extracted  a  fraction  — 

30%  —  of  the  remaining  oscillation  energy  at  each  half  cycle  of  the 
motion . 
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3.  CORRELATIONS  WITH  EXPERIMENTAL  DATA 

3.1  SUMMARY  OF  TEST  DATA  AND  COMPUTATION  RESULTS.  Tables  3.1.1 
thru  3.1.4  summarize  the  input  and  output  parameters  of  the  computa¬ 
tions  of  this  report.  These  tables  include  all  the  test  data  for  Spot 
and  White  Perch  obtained  on  the  1973  and  1975  test  series  in  the 
Chesapeake  Bay.  For  the  1973  test  data  (the  "500"  Shot  Numbers)  the 
fishes'  horizontal  ranges  and  depth  which  are  listed  in  the  tables 
are  the  preset  or  intended  values.1  On  the  1975  tests,  shots  783 
thru  786,  water  currents  caused  significant  deviations  in  horizontal 
range  and  depth.  For  these  shots  the  preset  values2  have  been 
corrected  as  described  in  Appendix  D.  (Tables  D-2  thru  D-7  list 
corrected  values  for  charge  depth  and  fish  and  pressure  gage  locations 
for  all  shots  of  the  1975  test  series.) 

The  ambient  pressure  (Tables  3.1.1  and  3.1.2)  is  the  total 
hydrostatic  head  at  fish  depth  calculated  for  fresh  water,  density  = 
1000  kg/m5,  and  nominal  atmospheric  pressure,  1.013  x  105  pascals. 

The  vapor  pressure  was  calculated  from  the  measured  water  temperature 
(1975  tests)  or  estimated  water  temperature  (1973  tests). 

The  peak  pressure  PMAX  in  pascals  and  time  constant  0  in 
milliseconds  are  calculated  from  the  similitude  relations  for 
pentolite 


PMAX 

0 


PCORR  x  5.62  x  10 7 

.1/3  /  R  \  0.23 


(•42) 


14 


0.084  WJ 


,W 


173, 


(3.1.1) 

(3.1.2) 


where  W  is  the  explosive  mass  in  kilograms,  R  is  the  slant  range  from 
the  charge  in  meters,  and  PCORR  is  an  empirical  constant  that  was  used 
to  adjust  Equation  3.1.1  for  batch-to-batch  variations  in  the  explo¬ 
sive  charges.  For  the  1973  test  series  PCORR  =  1.10  —  for  the 
1975  tests  PCORR  =  0.95. 


(Text  continued  on  nacre  32) 


TABLE  3.1.1  APPLIED  PRESSURE  PARAMETERS  AND  BLADDER  OSCILLATION  RESPONSE  FOR  SPOT 

(1  bar  ^  105  pascals  =  14.50  psi) 
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TABLE  3.1.2  APPLIED  PRESSURE  PARAMETERS  AND  BLADDER  OSCILLATION  RESPONSE  FOR  WHITE  PERCH 

(1  bar  =  10®  pascals  =  14.50  psi) 
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The  positive  duration  TPOS  (see  sketch  Figure  3.1.1)  was 
calculated  as  described  in  Appendix  C  for  the  1973  tests;  for  the 
1975  tests  measured  values  were  used.  The  negative  pressure  PNEG 
and  its  duration  DTNEG  were  measured  from  the  pressure  records. 

The  last  five  columns,  Tables  3.1.1  and  3.1.2,  list  calculated 
parameters  for  the  oscillatory  response  of  the  swim  bladder.  The 
number  of  positive  cycles  is  the  positive  duration  measured  in 
complete  cycles  of  oscillation  starting  from  the  initial  "at  rest" 
expansion.  The  number  of  negative  cycles  locates  the  final  return  to 
ambient  pressure  (end  of  negative  phase  —  see  Figure  3.1.1)  in  terms 
of  the  bladder  oscillation.  "Zero"  is  taken  as  the  last  expansion 
occurring  before  surface  cut-off.  See,  e.g.,  Figure  2.3.1.  Thus, 
unless  the  number  of  negative  cycles  is  greater  than  0.5,  no  final 
expansion  takes  place  during  the  negative  phase.  The  number  of 
negative  cycles,  (At/T)Np,  is  a  measure  of  the  contribution  of  the 
negative  phase  to  the  final  bladder  expansion.  The  relationship  is 
roughly  as  follows: 

(At/T)Np  £  0.5  No  contribution 

0.5  <  (At/T)Np  <  1.0  Partial  contribution 

1.0  <  (At/T)Np  Full  contribution 

The  reader  is  referred  to  Appendix  B  for  a  more  complete  description. 

The  bladder  oscillation  parameters  X,  Y,  Z  are  defined  as 

follows : 


X  =  -100  In  AMIN/A i  (3.1.3) 

Y  =  100  In  A MAX/A i  (3.1.4) 

Z  =  X  +  Y  =  100  In  A MAX/AMIN  (3.1.5) 
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where  In  is  the  natural  logarithm,  A-^  is  the  initial  at  rest  bladder 
radius,  and  AMIN  and  AMAX  are  the  smallest  and  largest  radii  during 
the  oscillatory  response.  In  the  cases  studied  so  far,  the  first 
compression  has  always  been  AMIN,  while  AMAX  generally  occurs  during 
or  following  the  negative  phase.  We  will  use  the  bladder  oscillation 
parameters  X,  Y,  and  Z  in  Sections  3.2  and  3.3  to  discuss  and  correlate 
the  observed  injuries  to  fish  specimens  on  our  tests. 

The  bladder  radius,  A^  —  Tables  3.1.3  and  3.1.4,  column  7  — 
is  the  radius  for  a  sphere  of  the  same  volume  as  the  fishes'  bladder 
at  fish  depth  for  a  fish  which  has  not  acclimated,  i.e.,  for  a  fish 
at  equilibrium  with  one  atmosphere  pressure  which  is  isothermally 
compressed  to  hydrostatic  pressure  at  fish  depth.*  Bladder  radius  A^ 
for  non-acclimated  fish  is  given  by 


A.  =  (A.) 

1  1  o 


^1.013  x  10  A  1/3 


) 


(3.1.6) ** 


where  p^  is  the  pressure  at  fish  depth  in  pascals  and  (A^)o,  the 
radius  at  one  atmosphere  pressure,  is  given  by 

( Af) q/l  =  0.033  (Spot)  (3.1.7) 


(Ai)o/L  =  0.055  (White  Perch)  (3.1.8) 


where  L  is  the  fork  length  of  the  fish.  The  constants  0.033  and 
0.055  were  selected  for  each  species  by  trial  and  error  so  as  to 
optimize  the  correlation  between  experimentally  observed  injuries  to 
the  test  specimens  and  the  calculated  bladder  oscillation  parameter 
Z  (e.g.,  Figure  3.2.1). 


*This  Is  how  it  was  done  on  the  explosion  tests.  It  would  have  been  prohibitively 
expensive  to  have  let  the  fish  acclimate  to  testing  depth,  because  of  the  time 
required.  This  problem  is  discussed  by  Gaspin,  Wiley,  Peters,  1976. 

**  For  the  calculations  of  this  report  the  exponent,  l/3y,  where  Y  =  1.4  —  for 
adiabatic  compression  —  was  inadvertently  used.  Since  the  computations  reported 
here  should  be  quite  insensitive  to  this  error,  they  were  not  redone. 
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The  number  of  fish  of  a  given  average  size  at  each  test  specimen 
location  and  summaries  of  the  injuries  evaluated  upon  post-shot 
disection  are  listed  in  the  last  5  columns  of  Tables  3.1.3  and  3.1.4 


Thus,  e.g.,  on  Shot  531  (Table  3.1.3)  at  the  76.2  meter  range 
position  there  were  10  Spot,  all  10  of  which  received  injuries  of 
level  1  or  greater,  7  of  which  received  injuries  of  level  2  or  greater 
5  of  which  received  injuries  of  level  3,  and  none  of  which  received 
level  4  injuries.  (There  were  no  injuries  of  level  5  observed  on 
either  the  1973  or  1975  tests.) 


Table  3.1.5  is  a  condensed  version  of  the  code  used  to  evaluate 
the  explosion  damage  to  the  fish  upon  post-shot  dissection.3 
Successive  injury  levels  are  of  increasing  severity  and  each  fish  is 
classified  to  fall  into  one  and  only  one  level  of  injury.  Figure 
3.1.2  depicts  this  fundamental  design  of  the  injury  code.  Thus,  the 
fundamental  physical  problem  of  this  report  is  to  estimate  test  condi¬ 
tions  for  the  transitions  from  one  injury  level  to  the  next.  We  will 
do  this  by  estimating  the  probability  of  events  such  as  "occurrence  of 
injury  level  3  or  greater"  as  a  function  of  the  calculated  damage 
parameter  Z  (Equation  3.1.5).  It  is  convenient  —  but  not  fundamental 
to  the  analysis  —  that  such  "cumulative"  injuries  are  also  what  is 
of  practical  interest  in  explosion  testing.  For  example,  the  event 
"injury  of  level  3  or  greater"  has  been  found  to  correlate  with  the 
observed  fish-kill*  on  underwater  explosion  tests.  And,  the  event 
"injury  of  level  2  or  greater"  has  been  considered  to  give  the  fish 
"little  chance  to  survive  predation."1 


3.  The  injury  evaluation  code  is  due  to:  Hubbs ,  C.  L. ,  Schultz,  E.  P.,  and 
Wisner,  R.,  Unpublished  preliminary  report  on  "Investigation  of  Effects  on  Caged 
Fishes  from  Underwater  Nitro-Carbo-Nitrate  Explosions,"  U.  of  California,  Scripps 
Institute  of  Oceanography,  1960.  The  complete  code  is  also  listed  in  reference  1 


*  The  term  "observed  fish-kill"  refers  to  the  dead  and  dying  fish  found  on  the 
surface  and  bottom  following  an  underwater  explosion. 
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TABLE  3.1.5 


CODE  FOR  EVALUATING  VISIBLE  FISH  DAMAGE 


Injury 

Level 


0  No  damage 


1  Light  hemorrhaging  in  tissues  covering  kidney 


2  Light  hemorrhaging  throughout  body  cavity,  some  kidney 

damage,  but  gas  bladder  intact 


3  Severe  hemorrhaging  throughout  body  cavity,  gross  kidney 

damage,  and  gas  bladder  burst 


4  Partial  break-thru  of  body  wall,  bleeding  about  anus 


5 


Ruptured  body  cavity,  internal  organs  scrambled  or  lost 
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However,  should  the  reader  require  for  example  the  probability, 

P (level  2  injury),  it  is  just  the  difference,  P (level  2  or  greater) - 
P  (level  3  or  greater). 

3.2  OBSERVED  INJURIES  TO  SPOT  VS  CALCULATED  OSCILLATION  PARAMETER 
Figure  3.2.1  shows  the  experimentally  observed  injuries  (as 
shown  in  each  row  of  Table  3.1.3)  to  Spot  on  the  1973  and  1975  tests 
plotted  as  a  function  of  Z  =  100  In  A MAX/ AMIN .  There  are  three  plots, 
one  for  each  of  the  injury  levels  1,  2,  and  3.  The  plotted  points 
represent  the  percent  of  fish  of  a  given  size  at  a  given  specimen 
location  receiving  injuries  of  the  indicated  level  or  greater. 

To  investigate  the  functional  dependence  of  the  observed 
injury  on  the  calculated  oscillation  parameter  Z  we  constructed 
Table  3.2.1  and  Figure  3.2.2.  Table  3.2.1  was  constructed  from 
Table  3.1.3  by  reordering  (sorting)  the  entries  in  order  of  increas¬ 
ing  value  of  Z  and  then  separating  the  new  table  into  groups  repre¬ 
senting  approximately  100  fish.  The  total  observed  injuries  for  each 
group  was  then  summarized  by  the  successive  row  entries  in  Table  3.2.1. 
Thus,  for  the  first  group  with  a  mid-range  value  for  2=54;  41%  of  the 
fish  received  injuries  of  level  1  or  greater;  11%,  level  2  or  greater; 
and  5%,  level  3  or  greater. 

Figure  3.2.2  shows  plots  of  the  averaqed  injury  data  tabulated 
in  Table  3.2.1.  The  solid  curves  are  drawn  by  eye  throuqh  the  data 
points.  The  experimental  data  for  level  3  injuries  to  Spot  exhibit 
less  scatter  than  that  for  levels  1  or  2 ;  and  also  less  than  for  all 
levels  of  injury  to  White  Perch  (Section  3.3).  Whether  or  not  this 
has  any  physical  significance  is  not  known. 

Also  shown,  as  a  dashed  line  in  Figure  3.2.2,  is  the  equation 

p  =  - - - l -  (3.2.1) 

1  +  EXP [-0.055 (Z-125) ] 

where  p  is  the  probability  of  observing  injury  of  level  3  or  greater. 
Equation  3.2.1  is  the  Cumulative  Logistic  Probability  Function  with 
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FIG.  3.2.1  OBSERVED  INJURIES  TO  SPOT  AS  A  FUNCTION  OF  Z 
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TABLE  3.2.1.  AVERAGED  OBSERVED  INJURY  PROBABILITY  FOR  SPOT  AS  A  FUNCTION  OF 

CALCULATED  OSCILLATION  PARAMETER  Z 


OBSERVED  CUMULATIVE  PROBABILITY 
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FIO.  3.2.2  OBSERVED  CUMULATIVE  INJURY  TO  SPOT  AS  A  FUNCTION  Of  / 
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t  lu>  parameter  values,  0.05f>  and  l-'1',  1  i  t  t  ed  to  the  experimental  data. 
Tho  paramotoi  \*v»  1  vu'  O.OS'l  was  old  a  i  nod  by  considei  inn  all  ot  the  Spot 
and  White  Porch  data,  levels  1,  and  1  (tom  the  two  test  set  ies. 

The  paramotoi  value  175  --  the  value  ot  7,  corrospondino  t  o  '>0* 
injury  probability  --  i s  the  maximum  likelihood  tit  value  to  t he 
level  1  injury  data  tot  Spot  . 

Kquat  ion  l . .  1  t  m  considered  It'  be  a  suit  ieiently  prooise 
representation  ot  the  data  tot  most  appl  ieat ions,  sue  It  as  those 
described  tit  Sect  ion  4  ot  this  t  epot  t  . 


t.l  OBSERVED  INJUR  1  I  S  TO  WHtTE  PERUI1  VS  OAU'Ul.ATED  OSO 1  1 .1  .AT ION 

PARAMETER 

F  inure  '.t.l  shows  t  ho  experimentally  observed  itt  juries  to 
White  Perch  on  t  lie  1*>7J  and  1*175  tests  (as  shown  in  each  row  ot 
Table  1.1.4)  p  let  t  let!  as  a  tunet  ion  ot  .  To  invest  ioate  the 
functional  dt'pendenc  ies  we  constructed  Table  t.i.1  and  K  inure  t .  ' . 
t  t  om  the  data  presented  in  Table  1.1.4.  The  plot  shown  in  Pin  tire 
1 .  t . for  level  l  injuries  exhibits  excessive  scatter.  Perhaps 
these  level  1  in  jut  ies  to  White  Perch  are  more  closely  related  tv' 
some  other  parametei  ot  the  bladdot  oscillation,  but  it  did  not 
appear  worthwhile  to  invest  inate  that  possibility. 

Fiuure  t.  t.  1  shows  the  avoranod  injury  data  tv'i  both  Spot  ami 
White  Perch  (Tables  1.7.1  and  t.  1.1)  plotted  as  a  tunet  ion  ot  ,  The 
curves  were  drawn  by  eye  throunh  t  hi'  combined  data.  These  plots 
support  the  viewpoint  ot  this  report  that  the  observed  in  jut  ies  t o 
both  species  of  fish  --  Spot  and  White  Perch  --  are  described  by  the 
same  daniane  functions.  (The  ooetfioionts  --  Equations  1,1.7  and 
1.1. H  --  for  the  effective  bladder  radii  used  in  the  dynamical 
response  calculations  are  difterent  tor  the  t wo  species  ot  tish, 
however . 1 
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FIG.  3.3.3  OBSERVED  CUMULATIVE  INJURY  AS  A  FUNCTION  OF  Z- 
SPOT  AND  WHITE  PERCH  COMBINED  DATA 
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3.4  COMPARISON  OF  SHOCKWAVE  IMPULSE  AND  CALCULATED  BLADDER 
OSCILLATION  AS  DAMAGE  PARAMETERS 

Yelverton  et  al1*  at  the  Lovelace  Foundation  working  in  a  test 
pond  with  eight  species  of  bladder  fish  at  depths  down  to  10  feet 
obtained  fish-kill  results  describable  by  an  Impulse  Damage  Parameter, 


Zx  =  80  In  [I/M1//3] 


(3.4.1) 


where  In  is  the  natural  logarithm,  I  is  the  impulse  of  the  positive 
portion  of  the  pressure  wave  in  psi-msec  and  M  is  the  mass  of  the 
fish  in  grams.  [The  arbitrary  constant,  80,  was  calculated  by  the 
author  to  correspond  to  the  value,  100,  in  the  Bladder  Oscillation 
Parameter  (Equation  3.1.5).] 

The  results  reported  by  Yelverton  et  al4  give  a  kill- 
probability 


P 


1 

1  +  EXP [-0.083 (Zj-132) ] 


(3.4.2) 


where  Z^.  is  the  Impulse  Parameter.  Equation  3.4.2  is  derived  from 
Figure  6  of  Yelverton  et  al  by  fitting  the  cumulative  logistic 
probability  distribution  to  their  results  for  observed  mortality  upon 
holding  the  fish  for  24  hours  after  the  explosion  test. 

Figure  3.4.1a  shows  the  level  3  injury  results  for  Spot  and 
White  Perch  (Tables  3.1.3  and  3.1.4)  plotted  as  a  function  of  the 
Impulse  Damage  Parameter.  Note  that  the  Impulse  Damage  Parameter 
does  not  describe  these  results . 

Figure  3.4.1b  shows  a  subset  of  the  level  3  injury  results 
for  Spot  and  White  Perch  plotted  as  a  function  of  the  Impulse  Damage 


4.  Yelverton,  J.  T.,  et  al,  ’’The  Relationship  between  Fish  Size  and  their  Response 
to  Underwater  Blast,"  DNA  Report  3677T,  Lovelace  Foundation,  1975. 
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Parameter.  In  Figure  3.4.1b  those  data  points  where  the  bladder 
response  calculation  showed  more  than  one  oscillation  during  the  posi¬ 
tive  phase  (Tables  3.1.3  and  3.1.4)  have  been  removed.  The  curve  is 
Equation  3.4.2,  the  kill-probability  result  from  Yelverton  et  ai 
(1975).  Thus,  their  24-hr  mortality  is  approximately  equivalent  to 
our  level  3  dissection  injuries;*  and,  perhaps  more  important,  the 
Impulse  Damage  Parameter  is  shown  to  describe  only  those  explosion 
geometries  where  the  charge  and/or  the  fish  are  at  shallow  enough 
depth  that  approximately  one  or  less  cycles  of  bladder  oscillation 
occur  before  surface  reflection  terminates  the  positive  pressure 
phase . 

Figure  3.4.1c  shows  for  comparison  all  of  the  level  3  injury 
results  plotted  as  a  function  of  the  Bladder  Oscillation  Parameter. 
The  curve  is  Equation  3.2.1  derived  from  the  Spot  level  3  injury 
data . 

3.5  DISCUSSION  OF  INJURY  CORRELATIONS 

It  is  not  accidental  that  both  the  bladder  oscillation 
parameter  and  the  impulse  damage  parameter  can  be  used  to  describe 
the  fish-kill  for  shallow  explosion  geometries.  To  see  how  this 
comes  about  we  examine  the  response  of  an  air  bubble  to  impulsive 
pressure  loading. 

"Let  the  pressure  be  applied  suddenly  and  let  it  disappear 
again  before  the  bubble  has  had  time  to  change  appreciably  in  size. 
Then,  the  bubble  will  begin  contracting  inward  at  a  certain  radial 
veloc i ty  v^  given  by 


v. 

l 


I 


PA, 


(3.5.1) 


*In  examining  dead  and  disabled  fish  which  have  been  collected  from  the  water 
surface  following  underwater  explosions  Martin  Wiley  and  Greig  Peters  (Chesapeake 
Biological  Laboratory,  Solomons,  Md.)  have  found  only  fish  having  injuries  of 
level  3  or  greater.  This  indicates  that  fish  having  received  lesser  injuries  do 
not  show  up  in  the  visible  fish  kill  although  they  mav  later  fall  victim  to 
predation . 
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where  p  is  the  density  of  the  water  and  I  = jp  dt,  the  applied 
impulse."5  By  the  methods  of  Appendix  A  we  can  now  calculate  the 
radius  ratio  AMAX/AMIN  for  the  oscillating  bubble.  See  Figure  3.5.1. 

The  solution  —  by  means  of  a  tabulated  function  —  is  of  the 

form 


AMAX 

AMIN 


FUNCTION 


(3.5.2) 


where  p^  is  the  initial  value  of  the  ambient  water  pressure. 


Since  yj  pV  does  not  vary  greatly  for  shallow  explosion 
geometries  and  since  in  most  bladder  fish  the  swim  bladder  comprises 
a  roughly  constant  fraction  of  the  total  volume  (about  6%),  Equation 
3.5.2  can  be  written 


AMAX 

AMIN 


FUNCTION 


(approximately) 


(3.5.3) 


where  we  have  substituted  M,  the  mass  of  the  fish,  for  the  volume, 
since  all  fish  are  approximately  neutrally  buoyant.  Kauation  3.5.3 
shows  that  the  Bladder  Oscillation  Parameter  and  the  Impulse  Damage 
Parameter  are  for  practical  purposes  equivalent,  for  the  special 
condition  of  shallow  fish  depth  and  impulsive  pressure  loading. 

Taken  together,  the  results  of  the  Lovelace  Foundation  for 
fish-kill  as  a  function  of  the  impulse  and  the  present  results 
described  in  terms  of  the  bladder  oscillation  parameter  (and,  in 
part  by  the  impulse  damage  parameter)  give  us  confidence  that  we 
have  achieved  a  correct  approximate  solution  to  the  problem  of 


5.  Kennard,  E.  H.,  "Radial  Motion  of  Water  Surrounding  a  Sphere  of  Gas  in 
Relation  to  Pressure  Waves,"  1943,  published  inVol.  II  of  "Underwater  Explosion 
Research,"  Office  of  Naval  Research,  1950. 
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predicting  fish-kill.*  The  non-ideal  appearance  of  the  injury 
correlations  shown  in  Figure  3.4.1b  and  c  are  no  doubt  due  to  the 
approximate  nature  of  these  two  damage  parameters.  That  the  impulse 
damage  parameter  can  be  used  to  describe  observed  injuries  for 
pressure  waves  of  more  than  1/2-cycle  in  duration  (Figure  3.4.1b)  is 
clear  evidence  of  its  approximate  nature.  Apparently,  for  shallow 
explosion  geometries  these  two  damage  parameters  approximate  some 
unknown  "true"  damage  parameter  about  equally  well. 

Additional  work  could  possibly  yield  a  better  damage  parameter 
than  either  the  Impulse  Damage  Parameter  or  the  Bladder  Oscillation 
Parameter.  However  such  additional  work  does  not  seem  justified 
because  at  present  our  ability  to  predict  fish-kill  is  limited  by 
other  more  critical  factors,  namely,  knowledge  of  the  pressure 
signature  in  the  presence  of  surface  and  bottom  effects  and  knowledge 
of  the  fish  species/size/density  distribution  in  the  vicinity  of  the 
explosion . 

Finally,  two  observations  are  noted: 

(1)  The  bladder  oscillation  solution  can  be  used  to  define 
the  useful  range  of  the  impulse  damage  parameter. 

(2)  For  shallow  explosion  geometries  fish-kill  computations 
using  the  simpler  impulse  damage  parameter  can  be  used  to  check 
computations  using  the  computed  bladder  oscillation  parameter. 


*We  are  also  confident  (based  on  unpublished  work  to  date)  that  the  Lovelace 
Foundation  results  (Yelverton  et  al,  1975)  can  also  be  described  in  terms  of 
the  bladder  oscillation  parameter.  To  do  so  is  not  a  routine  matter,  however, 
since  the  negative  phase  pressure  and  duration  at  each  fish  location  must  be 
calculated,  and  an  effective  bladder  radius  coefficient,  (Aj)0/L  (analogous  to 
Equations  3.1.7  and  3.1.8),  must  be  determined  by  trial-and-error  calculations 
for  each  fish. 
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4  APPLICATIONS 

4.1  FISH-KILL  CONTOURS 

Figure  4.1.1  shows  the  region  of  greater  than  50%  kill  pre¬ 
dicted  for  16-cm  long  Spot  in  the  vicinity  of  a  32  kilogram  pentolite 
explosion  at  9  meter  depth.  The  procedure  for  making  Figure  4.1.1  is 
as  follows.  First,  we  calculate  the  bladder  oscillation  parameter  Z 
for  an  array  of  horizontal  ranges,  x,  and  depths,  y;  i.e.,  Z  =  f(x,y) 
These  values,  Z  =  100  In (AMAX/AMIN) ,  are  tabulated  in  Table  4.1.1. 
Next,  assuming  that  injury  of  level  3  constitutes  kill,  we  read  from 
Figure  3.2.2  the  value,  Z,-q%  =  120,  for  50%  kill  probability.  Finally 
we  plot  in  Figure  4.1.1  the  regions  of  Table  4.1.1  having  Z  areater 
than  Z50%. 

To  calculate  Z  in  Table  4,1.1  we  first  calculated  an  approxi¬ 
mate  pressure-time  signature  by  the  method  described  in  Appendix  C. 
Then,  we  calculated  the  response  to  this  signature  as  described  in 
Section  2  and  in  Appendix  B,  except  that  these  computations  were  for 
undamped  oscillatory  bladder  response.* 

Figure  4.1.2  shows  predicted  regions  of  greater  than  10%,  50% 
and  90%  kill  for  21.5-cm  long  White  Perch  for  32  KG  pentolite  at  9 
meter  depth.  The  procedure  for  making  Figure  4.1.2  was  the  same  as 
described  for  Figure  4.1.1.  Values  of  Z1Q%  =  77  and  Zgo%  =  190  read 
from  Figure  3.2.2  were  used  to  obtain  the  10%  and  90%  kill  contours 
in  Figure  4.1.2.  in  Figure  4.1.2  contour  details  such  as  shown  in 
Figure  4.1.1  have  been  smoothed  out. 

Note  two  important  features  of  our  solution  to  the  fish-kill 
problem  which  are  evident  in  Figure  4.1.2: 

(1)  The  fish-kill  is  strongly  dependent  on  the  depth  of 

the  fish. 


*Due  to  "no  damping"  and  because  of  discrepancies  between  the  negative  pressures 
and  durations  calculated  by  Appendix  C  and  results  from  the  1975  test  series 
these  calculations  are  considered  approximate  —  but  quite  adequate  for  the 
purposes  of  this  report. 
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TABLE  4.1.1 

BLADDER  OSCILLATION  PARAMETER  Z  AS  A  FUNCTION  OF  FISHES'  DEPTH 

AND  HORIZONTAL  RANGE* 


Charqe:  32  KG  Pentolite  at  9  M  Depth 

Fish:  Spot,  18  cm  Fork  Lenqth 


HORIZONTAL  RANGE  (METERS) 


61 

91 

122 

152 

183 

213 

244 

0 

1.5 

29  9 

333 

131 

94 

52 

44 

39 

3.0 

302 

291 

139 

153 

111 

112 

72 

4.6 

310 

288 

111 

96 

131 

143 

98 

6.1 

228 

217 

194 

139 

102 

98 

122 

7.6 

244 

185 

193 

158 

- 

124 

128 

9.1 

174 

207 

150 

171 

131 

104 

119 

10.7 

228 

151 

135 

147 

137 

103 

98 

12.2 

156 

140 

149 

115 

136 

111 

89 

13.7 

205 

159 

140 

94 

109 

111 

90 

15.2 

149 

116 

95 

110 

92 

112 

99 

16.8 

188 

131 

102 

123 

75 

90 

100 

18.3 

142 

138 

132 

98 

90 

66 

87 

19.8 

175 

107 

91 

74 

103 

6  2 

67 

21.3 

136 

140 

86 

78 

90 

76 

55 

22.9 

163 

95 

106 

104 

68 

87 

60 

24.4 

133 

91 

92 

82 

59 

76 

65 

25.9 

110 

114 

80 

6  7 

71 

58 

75 

27.4 

128 

96 

95 

76 

85 

50 

66 

29.0 

125 

119 

82 

91 

57 

53 

50 

30.5 

118 

92 

74 

60 

54 

73 

42 

"These  values  for  Z  were  calctilated  using  approximate  predicted  pressure-time 
inputs  and  no  damping  of  the  calculated  bladder  motion. 
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(2)  Regions  of  less  than  50%  kill  —  since  they  encompass 
a  tremendous  volume  of  water  —  make  a  major  contribution  to  the 
total  number  of  fish  killed  by  an  underwater  explosion. 

These  two  features  are  examined  further  in  Sections  4.2  and  4.3. 

4.2  KILL  PROBABILITES  AS  A  FUNCTION  OF  DEPTH  FOR  DIFFERENT 
SIZED  FISH 

Figure  4.2.1  shows  the  variation  of  predicted  kill-probability 
with  fishes'  depth  at  a  fixed  horizontal  range  for  different  fish 
corresponding  to  three  different  sized  swim  bladders.  The  predicted 
kill-probability  is  obtained  from  calculated  Z-values,  e.g.. 

Table  4.1.1,  using  Equation  3.2.1  (dashed  line  in  Figure  3.2.2).  The 
approximate  length  of  Striped  Bass  corresponding  to  the  equivalent 
bladder  radii  of  the  bladder  oscillation  calculations  was  calculated 
from 


Bladder  Radius 
Fish  Length 


(4.2.1)* 


For  the  two  larger  sized  fish  in  Figure  4.2.1  the  maxima 
occurring  at  about  7  and  11  meters  depth,  respectively,  are  caused  by 
surface  cut-off  occurring  at  the  first  half-cycle  of  bladder 
oscillation.  At  shallower  depths  the  bladder  does  not  have  time  to 
respond  fully  to  the  positive  portion  of  the  explosion  wave;  while 
coincidence  of  the  first  half-cycle  and  surface  cut-off  amounts  to 
a  resonance  between  the  oscillatory  response  and  the  driving  pressure 
field  outside.  Thus,  at  shallow  depth  the  larger  fish  are  in  effect 
protected  from  harm  by  their  swim  bladders;  while  at  the  resonance 
depth  their  swim  bladders  "do  them  in." 


*The  value  .042  is  80%  of  that  experimentally  measured  from  a  sample  of  3 
Stripped  Bass.  The  value  80%  represents  a  guess  at  the  correction  to  go 
from  measured  radius  to  equivalent  computation  radius  for  this  fish. 
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4.3  TOTAL  FISH-KILL  AND  ITS  SPATIAL  DISTRIBUTION 

In  this  section  we  calculate  the  total  fish-kill  from  an 
underwater  explosion  and  its  distribution  as  a  function  of  horizontal 
range  and  depth,  assuming  a  uniform  spatial  distribution  of  a  single 
size  class  of  one  species  of  fish. 

Table  4.3.1  (top)  lists  the  bladder  oscillation  parameter  Z 
for  21.5  cm  White  Perch  calculated  over  a  cylindrical  grid  surround¬ 
ing  a  32  kilogram  charge  fired  at  9  meter  depth.  Table  4.3.1 
(bottom)  lists  the  corresponding  kill  probabilities  from  Equation 
3.2.1. 

To  calculate  the  fish-kill  we  must  know  which  fish  are 
where;  and,  in  general,  this  is  impossible.  Nevertheless,  to  gain 
insight  as  to  the  potential  of  a  given  explosion  geometry  for  killing 
fish  it  is  useful  to  assume  a  fictitious  fish  density  distribution. 

For  this  purpose  we  will  assume  a  nominal  uniform  density  distribution 
of  1  fish  per  thousand  cubic  meters  of  water. 

To  calculate  the  fish-kill  we  must  integrate  the  product, 
fish  density  x  kill  probability,  over  the  region  surrounding  tne 
explosion.  Thus,  for  this  example  we  compute  the  products  of  "fish 
present"  and  "kill  probability"  (Tables  4.3.2  and  4.3.1)  to  get  the 
"fish  killed"  in  each  mesh  volume,  Table  4.3.2  (bottom)  —  each  mesh 
consisting  of  a  cylindrical  annulus  100  feet  thick  and  5  feet  in  depth. 
Summing  the  distribution  of  "fish  killed"  shown  at  the  bottom  of 
Table  4.3.2  results  in  a  total  kill  of  1245  fish  out  of  6228  fish 
present.  Note  that  we  have  neglected  200  fish  present  in  the  150-ft 
radius  cylinder  of  water  containing  the  charge  as  well  as  all  fish 
present  beyond  a  horizontal  range  of  850  feet.  Nevertheless,  we 
have  described  the  most  significant  portion  of  the  distribution  of 
fish  killed  by  this  underwater  explosion. 

It  is  apparent,  for  example  that  if  fish  were  present  only 
at  depths  greater  than  70  or  80  feet  relatively  few  fish  would  be 
killed.  Likewise,  if  fish  (of  this  particular  size)  were  present 
only  within  some  10  feet  of  the  surface  the  kill  would  be  relatively 
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TABLE  4.3.1  BLADDER  OSCILLATION  PARAMETER  AND  KILL  PROBABILITY  AS  A 
FUNCTION  OF  FISHES'  DEPTH  AND  HORIZONTAL  RANGE 


Charge:  32  KG  Pentolite  at  9  M  Depth 

Fish:  White  Perch,  21.5  cm  Fork  Length 


Horizontal  Range  (feet) 


200 

300 

400 

500 

600 

700 

800 

5 

- 

251 

208 

41 

41 

35 

10 

7 

10 

- 

220 

21  1 

119 

73 

32 

30 

28 

IS 

- 

206 

215 

156 

100 

62 

31 

26 

20 

- 

287 

197 

131 

129 

86 

79 

58 

25 

206 

180 

147 

114 

109 

74 

69 

TO 

- 

160 

175 

1  30 

146 

94 

93 

64 

35 

- 

172 

187 

117 

120 

113 

88 

81 

Fish 

40 

- 

211 

163 

120 

97 

102 

97 

96 

BLADDER 

45 

- 

168 

125 

135 

89 

94 

99 

84 

Depth 

50 

- 

144 

110 

139 

94 

76 

82 

87 

OSCILLATION 

(feet) 

55 

- 

133 

111 

121 

99 

70 

77 

73 

PARAMETER 

60 

- 

177 

133 

91 

1  04 

76 

62 

68 

65 

- 

139 

141 

83 

1  08 

72 

58 

55 

70 

- 

132 

110 

81 

94 

86 

63 

51 

75 

- 

128 

98 

79 

70 

91 

60 

49 

“0 

- 

154 

06 

98 

67 

87 

66 

54 

85 

- 

107 

1  10 

105 

63 

76 

70 

52 

90 

- 

120 

119 

81 

61 

56 

75 

57 

95 

- 

125 

76 

66 

69 

54 

65 

61 

100 

122 

08 

71 

84 

51 

56 

65 

200 

300 

400 

500 

600 

700 

800 

5 

- 

1.00 

.99 

.01 

.01 

.01 

.00 

.00 

10 

- 

.99 

.99 

.42 

.05 

.01 

.01 

.00 

15 

- 

.99 

.99 

.85 

.20 

.03 

.01 

.00 

20 

- 

1.00 

.98 

.58 

.55 

.10 

.07 

.02 

25 

- 

.99 

.95 

.77 

.35 

.29 

.06 

.04 

l 

30 

- 

.91 

.94 

.57 

.76 

.15 

.15 

.03 

35 

- 

.93 

.97 

.39 

.43 

.34 

.12 

.08 

Fish 

40 

- 

.99 

.89 

.43 

.18 

.22 

.18 

.17 

45 

- 

.91 

.50 

.63 

.12 

.15 

.19 

.09 

KILL 

Depth 

50 

- 

.74 

.30 

.68 

.15 

.06 

.09 

.11 

55 

- 

.61 

.32 

.45 

.19 

.05 

.07 

.05 

PROBABILITY 

,  I  e  e  ) 

60 

- 

.95 

.61 

.13 

.24 

.06 

.03 

.04 

65 

- 

.68 

.71 

.09 

.28 

.05 

.02 

.02 

70 

- 

.60 

.30 

.08 

.15 

.10 

.03 

.02 

1  I 

75 

- 

.54 

.18 

.07 

.05 

.13 

.03 

.02 

80 

- 

.83 

.10 

.18 

.04 

.11 

.04 

.02 

85 

- 

.27 

.30 

.25 

.03 

.06 

.05 

.02 

90 

- 

.43 

.42 

.06 

.03 

.02 

.06 

.02 

9 

95 

- 

.50 

.06 

.04 

.04 

.02 

.04 

.03 

100 

.46 

.12 

.05 

.09 

.02 

.02 

.04 
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TABLE  4.3.2  FISH  PRESENT  AND  NUMBER  KILLED  AS  A  FUNCTION  OF  WATER  DEPTH 
AND  HORIZONTAL  RANGE 

Charge:  3 2  KG  Pentolite  at  9  M  Depth 

Pish:  White  Perch,  21.5  cm  Fork  Length 

Fish  Density  ^ 

Distribution:  Uniform  Density  =  1  x  10-:S  Fish/Meter'5 
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light.  However,  for  fish  of  this  size  between  depths  of  20  and  60 
feet  the  kill  will  be  extensive  out  to  a  range  of  about  800  feet; 
and  if  a  large  school  several  hundred  feet  in  diameter  with  a  density 
say,  of  one  fish  per  cubic  meter  were  within  this  region,  a 
disastrous  fish-kill  would  result.  (Such  an  occurrence  is  unlikely, 
however,  as  it  is  Navy  policy  to  delay  operations  if  a  school  of 
fish  is  observed  in  the  vicinity  of  the  test  site.) 
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APPENDIX  A 

RESPONSE  OF  GAS  BLADDER  TO  UNDERWATER  EXPLOSION  PRESSURES: 
SQUARE  STEP  APPROXIMATION6 

ABSTRACT 

This  appendix  presents  the  theory  of  the  radial  motion  of 
a  gas-filled  spherical  bubble  responding  to  step  changes  in  the  water 
pressure . 


The  theory  is  used  to  estimate  the  motion  of  a  fish-bladder 
sized  bubble  subjected  to  the  pressure  signature  from  an  underwater 
explosion.  Some  important  qualitative  results  as  well  as  some  crude 
quantitative  results  are  readily  obtained.  For  example: 

(1)  For  a  simple  positive  step  which  returns  to  ambient 
pressure,  the  maximum  bladder  extension  depends  on  both  the 
magnitude  and  duration  of  the  pressure  pulse.  If  cut-off  occurs 
when  the  bladder  (undamped  oscillation)  is  at  maximum  size,  no 
extension  beyond  its  initial  at-rest  size  occurs.  If  cut-off  occurs 
at  any  other  time,  some  over-extension  of  the  bladder  occurs. 

(2)  A  sudden  decrease  in  the  ambient  pressure  will  always 
result  in  some  over-extension  of  the  bladder.  If  this  negative 
pressure  (below  ambient)  follows  a  positive  step  with  cut-off 
occurring  at  maximum  size,  the  over-extension  is  due  solely  to  the 
negative  portion  of  the  pressure  signature;  if  not,  the  over-extension 
is  due  to  both  the  positive  step  and  the  negative  pressure  following 
cut-off.  The  effects  are  additive  in  a  non-linear  manner  described 

in  this  appendix. 


6.  Appendix  A  was  formerly  given  internal  distribution  as  NOLTN  10205  ,  Goertner  , 
1974,  "Response  of  Air  Bubbles  to  Underwater  Explosion  Pressures:  Square  Step 
Approximation  with  Application  to  Fish-bl  aider-cavity  Response". 
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I.  ON  THE  MOTION  OF  AN  AIR  BUBBLE  SUBJECTED  TO  A  SERIES  OF  STEP 
CHANGES  IN  THE  OUTSIDE  WATER  PRESSURE 


1.  Consider  a  spherical  bubble  at  rest  in  water  at  the 
initial  value  of  the  ambient  pressure.  Let  A-^  be  its  initial 
radius  when  at  rest.  Its  internal  air  pressure  is  the  initial 
ambient  pressure  p^.  At  t  =  o,  the  ambient  pressure  jumps  to  a  new 
value  p^ ;  and  the  bubble  is  then  oscillating.  The  situation  for  the 
air  pressure  inside  the  bubble,  trying  to  adjust  to  the  new  outside 
pressure,  is  sketched  in  Figure  A-l. 


2.  In  this  Appendix  we  concern  ourselves  with  the  motion 
of  a  bubble  subjected  to  a  succession  of  jumps  to  new  ambient 
pressures  pn.  At  each  ambient  pressure  the  bubble  motion  is 
described  by  the  following  differential  equation J 


|  (i^/OA3)  A2  +^P0A3  +  E  (A)  =  Y 


Kinetic  Energy  Potential 

of  Energy  of 

Surrounding  Water  Pushed-back 

Water 


Total  Energy 


Internal  Energy 
of  Air 

Inside  Bubble 


7.  Cole,  R.  H.  ,  19U8,  "Underwater  Explosions",  Princeton  University  Press, 
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Pi  =  Initial  outside  water  pressure 
p^  =  New  value  of  outside  water  pressure 
P^  =  Bubble  air  pressure  at  maximum  size 
Pq  =  Bubble  air  pressure  at  "Equilibrium" 

=  Pi 

P  =  Bubble  air  pressure  at  minimum  size 
m 

FIG.  A-1  BUBBLE  AIR  PRESSURE  AND  OUTSIDE  WATER  PRESSURE  AS  A 
FUNCTION  OF  TIME  (SKETCH) 
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whoro 


P (A)  x  V(A) 


Bubble  Radius 

,  the  radial  velocity 
Aii-  Pressure  Inside  the  Hubble 


Hubble  Volume 


4  n  .3 


Kxponont  in  Adiabatic  P-V  Relationship  1.4 
(for  air) 

Water  Density  constant 

Ambient  Water  Pressure  constant 

Total  Kneruy  ot  Hubble  Oscillation  constant 


3.  The  oscillatory  motion  described  by  liquation  Al  has 
much  in  common  with  the  osclllat ion  ot  a  simple  mass-sprinq  system 
the  difference  is  that  for  the  bubble  oscillation  the  "mass"  chanues 
with  time  and  t  lio  "spring"  is  non-linear.  The  overall  behavior  of 
oscillating  gas  bubbles  was  worked  out  in  considerable  detail  by 
Snay  and  Christian. 

They  rewrote  liquation  Al  in  the  non-dimensional  form 

a  1  A  ‘  t  a  ’  +  k  a  ^ 1  *  ^  ”  1 )  « 1  ( A  1 ) 

where  a  is  the  non-dimensional  radius  a  A/1,  and  a  is  the  derivative 
of  "a"  with  respect  to  the  non-d imens i on a  1  t  ime  t  '  t/C.  The 
parameters  k  and  >  characteristic  the  bubble  mid  ion.  "k"  can  be 
expressed  by 


KW-V) 


p„  h  n 


whore  P  is  the  internal  pressure  at  either  extremum  ot  the  bubble 
oscillat  ion  M  refers  to  the  maximum  bubble  volume  (P.,  minimum 

N 

pressure),  m  refers  to  the  minimvtm  bubble  volume  (P^  maximum 
pressure)  —  and  P0  is  the  ambient  hydrostatic  pressure. 

fi  Snsy i  8i  0. t  ami  S>  a.  ChriitiiBi  1958 •  utadimtw  Explosion  PImdomm: 

Tlu-  Psirnmet  rv;i  of  n  Nen-MI  ^r«t  I  uk  Rnl'Mr  One  i  1  1  a t  I n^  in  an  InoompreimlMe 
Mr'.t  i  um"  ,  MAVOHh  Hr port  ,'l*  Vf 
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4.  The  lenqth  scale  factor  L  and  the  time  scale  factor  C 
are  defined  by 


(A5 ) 

(A6) 


where  Y  is  the  total  energy  (Equation  Al)  and  p  is  the  density  of 
the  water. 


5.  Snay  and  Christian  in  their  report  have  determined  the 
overall  behavior  of  the  oscillatory  solutions  to  Equation  A3  for  y 
in  the  range,  1.0  <  y  <1.5.  Table  Al  (below)  list  some  of  their 
results  for  y  =  1.4.  The  dimensionless  quantities  —  a^  a^  and  T' 
—  are  the  maximum  radius,  the  minimum  radius  and  the  period  of 
oscillation  of  the  bubble,  respectively. 


TABLE  A-l 

Bubble  Parameters  for  y  =  1.4 


k 

aM 

am/aM 

T  ' 

T  '  Calculated 

From  Equation  A32 

0.0 

1.000 

0.0 

1.492 

1.476 

0.05 

0.982 

0.085 

1.505 

1.496 

0.10 

0.963 

0.152 

1.522 

1.516 

0.15 

0.942 

0.220 

1.538 

1.536 

0.20 

0.920 

0.288 

1.557 

1.556 

0.25 

0.893 

0.364 

1.577 

1.576 

0.30 

0.863 

0.447 

1.596 

1.596 

0.35 

0.823 

0.548 

1.616 

1.616 

0.40 

0.768 

0.699 

1.636 

1.636 

0.432 

0.659 

1.000 

1.650 

1.649 
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6.  In  what  follows,  we  will  use  both  Equations  A1  and  A3. 
Transitions  from  one  outside  pressure  level,  p^,  to  another,  p  ^ , 
are  simply  described  using  Equation  Al.  And,  Snay  and  Christian's 
results,  e.g.,  those  listed  in  Table  I,  which  were  computed  from 
Equation  A3,  are  convenient  for  computing  the  bubble  behavior  within 
a  given  state. 

7.  Boundary  Conditions  at  Pressure  Jumps.  We  now  consider 


a  step  change  in  outside  pressure  from  p^  to  Pn+1' 


Let  V  be  the 
c 

n 


bubble  volume  at  the  instant  of  change,  and  y  and  y  , .  be  the  old 

th  n  n+1 

and  new  total  energies.  For  the  nr  state  we  rewrite  Eauation  Al  as 


(4 tt  p  A3  \a2  +  4up  A3  + 

-  )  -p" 


E (A)  =Y 


n 


(A7) 


Now,  since  the  oscillating  system  described  by  (A7)  has  finite  mass, 
neither  A  nor  A  can  change  impulsively,  i.e., 


AA  =  0 
AA  =  0 


(A8) 

(A9) 


at  each  pressure  jump.  Thus,  the  first  and  third  terms  of  (A7)  do 
not  change  at  the  jump,  and  the  change  in  the  total  enerqy  is  given 
by 


n+1 


=  Yn  +  -  Pj  V 


n+1 


(A10) 


where  V  is  the  bubble  volume  at  the  time  of  the  jump, 
n 

8.  Equations  (A7) ,  (A8),  (A10)  completely  specify  the 

motion  in  the  nfc^  state  given  suitable  initial  conditions.  In  the 
present  case  the  initial  or  zeroth  state  is  specified  by 

(All) 

(A12) 


A  =  constant  =  A^ 


P(A^)  =  constant  =  p^ 
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In  particular,  the  total  energy  and  the  outside  pressure  pn  are 
sufficient  to  specify  the  amplitude  of  the  oscillation.  And,  if  the 
radius  function  A(t)  is  known,  the  phase  will  then  be  determined  by 
the  jump  conditions,  AA  =  0  and  AA  =  0. 

9.  Some  General  Equations  for  the  n*'*1  State.  We  now  list 
some  further  equations  that  can  be  used  to  calculate  the  parameters 
of  the  bubble  oscillation  in  the  n*"*1  state. 


(A13) 


(A14) 


(A15) 


(A16) 


AMIX 

AMAX 

n 


P  (A) 


A 


n 


T 


n 


Tabulated  Function  (k  ) 


t(kn)  *  C„ 


( A17 ) 

(A18) 

(A19 ) 

( A20) 
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Equation  A19  gives  the  equilibrium  radius  A  (radius  after  the 

t  hi 

oscillation  has  damped  out)  for  the  bubble  in  the  nu  state.  Equation 
A19  is  derived  from  (A18)  by  setting  P(A)  equal  to  p^ . 

10.  Bubble  Motion  in  State  1.  Consider  the  step  change 
in  outside  pressure  from  p^  to  p^.  In  accordance  with  the  jump 
conditions  (A8)  and  (A9)  and  the  initial  conditions  (All)  and  (A12) 
we  set  A  =  0  and  A  =  A^  in  (A7)  to  obtain  the  total  energy 


Y1  *  <Pl  +  Vi  (A21) 

where  we  have  also  used  Equation  A2  to  evaluate  the  internal  energy 
E(Ai)  . 

11.  Since  A  =  0,  the  oscillation  in  state  1  starts  at  an 

extremum.  Whether  this  extremum  is  a  maximum  volume  V_,  or  minimum 

M 

is  determined  by  p^  being  "greater  than"  or  "less  than"  p^. 

If  p^  >  p^ ,  then 

AM2  =  Ai  (A22) 

\  =  Pi  (A23) 

IF  p^  <  p^,  then 

A  =  A.  (A24) 

m,  l 


(A25) 
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12.  Since  the  oscillation  in  state  1  starts  at  an  extremum, 
the  bubble  parameter  k  can  easily  be  computed  from  Eauation  A4.  From 
(A4)  together  with  (A23)  and  A25  we  get 


14.  Calculation  of  the  Bubble  Period.  The  bubble  period 
of  oscillation  T  is  the  product  of  the  dimensionless  period, 
t  =  F(y,k),  given  in  Table  I  and  the  time  scale  factor  C  (Equations 
A15  and  A14) 


=  T  '  x  C  (A29) 

Combining  (A14 ,  (A21)  and  (A27) ,  we  write  for  the  length  scale  factor 

L1  =  Ai  (1  +  xi)3/3  (A30) 
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And,  combining  (A15) ,  (A27)  and  (A30) ,  we  write  the  time  scale 

factor  as 


V3 

(1  +  XL)  A. 


( A31 ) 


where  p  is  the  water  density. 

15.  For  our  purposes  the  dimensionless  period  T1,  is 
conveniently  approximated  by  the  linear  function 

T  '  =  1.476  +  0.4  k  (A32) 

Values  of  T  '  calculated  by  Equation  A32  are  listed  in  Table  1 
alongside  the  tabulated  function.  The  approximation,  Equation  A32, 
is  good  to  better  than  0.5%  for  values  of  k  >  0.1. 

16.  Finally,  combining  (A26) ,  (A27) ,  (A29),  (A31)  and  (A32) 

we  write  for  the  period  of  oscillation 


T 


1 


(1  +  Xj_) 


1/3 


(1.476  +  0.4k:)  A i 


(A33) 


The  physical  quantities  in  Equation  A33  may  be  expressed  in  any 
consistent  set  of  units.  For  English  Engineering  Units  we  have 

T^  =  Bubble  Period  (SECONDS) 

p  =  Density  of  the  Water  (SLUGS/FT3) 

.  .  2 
_  Initial  Value  of  Ambient  Water  Pressure  (LBS/FT  ) 

"i  [pressure  in  LBS/FT^  =  (pressure  in  PSI)  x  144) 

A^  =  Initial  Radius  of  Bubble  (FT) 


A-10 
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17.  Simplifying  Assumption.  To  proceed  further  with  the 

general  case  we  would  now  need  the  bubble  radius  as  a  function  of 

time  in  order  to  apply  the  jump  equations  (A8) ,  (A9)  and  (A10)  at 

the  instant  the  outside  pressure  changes  to  its  new  value,  p2 . 

Although  not  particularly  difficult,  such  computations  are  beyond 

the  limited  scope  of  this  study.  Consequently,  at  this  point  we 

introduce  the  simplifying  assumption  that  all  pressure  changes  occur 

at  an  instant  when  A  =  0,  i.e.,  either  at  an  extremum  or  at 

equilibrium  bubble  radius  after  the  oscillation  has  completely 

damped  out.  Pressure  jumps,  p^  to  pn+^,  which  happen  to  occur  at 

bubble  extrema  give  the  greatest  and  also  the  least  amplitudes  of 

t  hi 

oscillation  in  the  n+l11  state  (V  =  V,,  in  Equation  A10)  .  Thus, 

c  Mm 

even  when  the  pressure  jumps  occur  at  other  times  than  the  extrema, 
we  can  place  exact  upper  and  lower  limits  on  the  amplitude  of 
oscillation  in  each  state. 

18.  Bubble  Motion  in  State  2.  In  accordance  with  our 

simplifying  assumption  we  will  assume  the  pressure  jump,  p^  to  p2 , 

occurs  at  a  time  t  such  that 

C1 

tc  -  N  x  j  Tj_  (A34 ) 

where  N  is  any  positive  inteqer.  Alternatively,  t  may  be  taken 

C1 

large  enough  that  the  bubble  oscillation  has  damped  out,  i.e., 

—  l/3y 

Ac  ~  =  (Pj/P^)  the  equilibrium  radius. 

19.  Let  V  ,  A  ,  P  be  the  volume,  radius  and  air 

C1  C1  C1 

pressure,  respectively,  of  the  bubble  at  time,  t  .  If  the  new  out- 

C1 

side  pressure  p„  is  higher  than  P  ,  the  air  pressure  inside  the 

C1 

bubble,  the  new  oscillation  begins  at  maximum  volume;  and  if  p2  is 

lower  than  P  ,  the  oscillation  begins  at  minimum  volume.  Or, 

C1 

restated  we  have 
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21.  Bubble  Motion  in  State  n.  The  corresponding  equations 
for  the  n1"  state  can  easily  be  obtained  from  those  for  state  2 
by  substituting  "n's"  for  "2's"  and  "n-l's"  for  "l's".  The  only 
apparent  exception  is  that  the  equivalent  equation  for  (A34)  should 
read 

T  =  T  +  N  x  i-T  (A42) 

c  c  .  2  n 
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II.  ON  THE  MOTION  OF  A  SWIM-BLADDER-SIZE  AIR  BUBBLE  SUBJECTED 
TO  AN  EXPLOSION  GENERATED  PRESSURE  WAVE 

22.  At  the  present  time  we  ask  such  questions  as:  "Can 
a  nearby  explosion  kill  or  seriously  damage  fish  by  extending  their 
swim  bladders?"  Prodded  by  this,  let  us  calculate  the  behavior  of 

a  bladder-sized  air  bubble  subjected  to  an  explosion  pressure 
signature — one  which  kills  many,  but  not  all,  of  the  fish  present. 

The  pressure  signature  shown  in  Figure  A-2  is  derived  from  a  experi¬ 
mental  test  condition  (Shot  525,  Chesapeake  Bay  Tests,  July-Aug  73) 
and  appears  to  meet  this  criterion.  It  corresponds  to  a  fish  cage 
location  where,  upon  subsequent  dissection,  all  ten  of  ten  5.4"  long 
Spot  and  six  of  ten  7.2"  long  White  Perch  were  judged  to  have 
received  lethal  damage.  For  this  calculation  we  will  further 
approximate  the  outside  pressure  signature  with  the  square-stepped 
one  which  is  also  shown  in  Figure  A-2. 

23.  To  get  a  better  feel  for  physical  phenomena,  we  will 
do  the  calculation  in  three  stages.  First,  we  will  calculate  the 
response  to  the  positive  portion,  a  square-step  which  returns  to  the 
initial  ambient  pressure  level;  then,  the  response  to  a  simple  step 
decrease  in  outside  pressure  to  the  level  of  the  explosion  induced 
underpressure;  and  finally,  the  combined  response  to  the  square- 
stepped  signature  sketched  in  Figure  A-2. 

24.  We  will  take  our  bladder-sized  air  bubble  from  a  5.4" 
long  Spot,  and  start  the  calculation  with  an  initial  (at  rest)  bubble 
radius  given  by* 


Aj[  =  0.0395  L 
Footnote  on  page  A-16 


(A43) 
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where  L  is  the  fork  length  of  the  fish.  For  our  5.4"  Spot  this  gives 
an  initial  bubble  radius 


Ai  =  0.0178  ft 


(A44) 


The  initial  ambient  pressure  p^  =  32  osi,  which  corresponds  to  a  depth 
of  40  feet  in  fresh  water.  The  magnitude  of  the  pressure  step, 

Pi  “  is  50  psi;  and  its  time  duration  TPOS  =2.3  milliseconds. 

The  period  T^  at  p^  =  82  psi  turns  out  to  be  0.569  milliseconds,  so 
that  the  oscillation  in  State  1,  which  starts  at  maximum  bubble  size 
am  =  g°es  through  2.3/0.569  =  4.04  periods  of  oscillation  by 

the  time  cut  off  occurs. 


(Footnote  from  page  A-14) 

*This  gives  a  rough  approximation  to  the  bladder  volume  of  a  5.4"  lone 
fish  at  40ft-depth  and  500  feet  distant  from  a  68-lb  pentolite  charge 
(Shot  525) .  Equation  A43  was  conjured  up  by  the  author  for  a 

t  h 

"nominal  fish".  It  represents  the  radius  of  a  sphere  of  1/20  of 

g 

the  estimated  volume  of  the  fish  (after  Weston) .  The  volume  of  the 
fish  is  taken  as  an  ellipsoid  whose  semi-major  axes  (A,B,C)  bear  the 
relation,  A  =  0.38  L,  B  =  0.30  A,  C  =  0.25  B,  where  L  is  the  overall 
length  of  the  fish.  The  above  coefficients  are  average  values  taken 
from  crude  measurements  of  sketches  of  fish  shown  by  Lippson. 

9.  Weston,  D.  E.  ,  1966,  "Sound  Propagation  in  the  Presence  of  Bladder  Fish", 
published  in  Vol.  II  of  "Underwater  Acoustics",  edited  by  V.  M.  Albers, 

1967,  Plenum  Press 

10.  Lippson,  A.  J. ,  1973,  "The  Chesapeake  Bay  in  Maryland  —  an  Atlas  of 
Natural  Resources",  John  Hopkins  University  Press 
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25.  At  cut-off  time,  TPOS,  the  measured  outside  water 

pressure  drops  to  about  13  psi  (Figure  A2)  .  Table  A2  lists  limiting 

case  results  for  the  final  maximum  radius  (in  State  2)  which  were 

calculated  from  the  above  parameter  values  by  assuming  cut-off  times 

TPOS  occurinq  at  A  ,  A  ,  and  also  sufficiently  late  that  the 

Mm  J 

oscillation  has  damped-out  (so  that  A  =  constant  =  A) . 

26.  Two  sets  of  values  for  A^  /A^  are  given  in  Table  A2 , 

one  for  a  final  outside  pressure  P2  =  pi  =  32  psi,  the  other  for 

the  measured  value  of  =  13  psi.  The  results  listed  in  the  top 

row  of  Table  A2  (cut-off  at  bubble  maximum)  are  the  same  as  if  the 

positive  outside  overpressure,  p^  -  ,  had  never  occurred.  In 

other  words,  the  upper  right  hand  value  of  A  /A.  also  applies  to  a 

M2  1 

simple  step  decrease  in  the  ambient  pressure  p.  to  a  new  value, 

P2  =  13  psi. 


27.  It  is  apparent  from  these  crude  results  that  both 
positive  and  negative  excursions  of  the  outside  pressure  significantly 
influence  the  final  overshoot;  and  that  time-of-occurrence  of  cut-off 
relative  to  the  phase  of  the  bubble  oscillation  cycle  is  of  critical 
importance  for  positive  excursions.  It  also  looks  like  a  good  bet 
that  underwater  explosions  do  kill  or  cause  serious  damage  to  nearby 


fish  by  extending  their  swim  bladders. 


28.  Figures  A3  and  A4  sketch  the  Qualitative  behavior  of 
the  inside  air  pressure  and  the  bubble  radius  as  functions  of  time 
for  the  six  cases  listed  in  Table  A2 . 


TABLE  A-2  DEPENDENCE  OF  RADIUS  OVERSHOOT  A^  /Aj  ON  CUT-OFF 
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PRESSURE  VERSUS  TIME 


RADIUS  VERSUS  TIME 


PRESSURE  VERSUS  TIME 


FIG.  A-4  AIR  BUBBLE  RESPONSE  TO  A  SQUARE  WAVE  FOLLOWED  BY  AN  UNDERPRESSURE 
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III.  DISCUSSION 


29.  After  the  above  analysis  was  completed,  George  Young 
brought  to  my  attention  a  similar  analysis  done  in  1943  by 
E.  H.  Kennard  (in  connection  with  the  proposed  use  of  screens  of 

5 

bubbles  as  a  protective  device  against  explosion  shock  waves).  His 
development  and  presentation  of  the  incompressible  bubble  theorv 
(pp.  190-191,  Ref.  5)  is  verv  convenient  for  calculating  bubble 
motions  such  as  the  example  used  in  this  note,  i.e.,  response  to 
jumps  in  constant  outside  pressure  states  restricted  to  those  instants 
when  the  radial  velocity,  A  =  0.  I  redid  the  calculations  summarized 
in  Table  A2  using  Kennard' s  "bubble  theory"  and  the  bubble  radii  and 
periods  of  oscillation  were  in  all  cases  within  1%  of  my  previous 
values  calculated  using  the  equations  presented  in  this  Appendix. 


A-21 


NSWC/WOL  TR  76-155 


IV.  CONCLUSIONS 


30.  I  hope  that  the  analysis  presented  here  (or  some 
further  development  of  it)  will  prove  useful  in  correlating  observed 
fish  injury  and  kill  with  incident  explosion  pressure  signatures. 

The  present  analysis  is  but  a  crude  first  cut,  however;  and  if 
this  approach  to  understanding  fish  damage  should  appear  fruitful, 
considerable  further  development  will  probably  be  desirable  and 
perhaps  essential  to  achieving  a  useful  correlation.  Nevertheless, 
some  things  are  already  clear.  For  example: 

(1)  By  itself,  a  step  increase  in  the  ambient  pressure 
cannot  extend  the  bubble/bladder  beyond  its  original  size,  since  the 
bubble  oscillates  between  its  original  volume  and  some  minimum  volume 
which  depends  on  the  magnitude  of  the  incident  shock  pressure. 
Extensions  beyond  the  initial  at-rest  volume  can  only  occur  after 
return  to  ambient  pressure  or  below. 

(2)  For  a  simple  positive  step  pressure  which  returns  to 
the  ambient  level,  the  maximum  bubble/bladder  extension  depends  on 
both  the  magnitude  and  duration  of  the  pressure  pulse.  If  cut-off 
occurs  at  the  instant  of  maximum  bubble  volume,  the  bubble  is 
returned  to  its  initial  at-rest  state,  having  undergone  no  extension 
beyond  this  initial  at-rest  state.  If  cut-off  occurs  at  any  other 
time,  some  extension  beyond  the  initial  at-rest  volume  occurs.  For 
a  given  incident  pressure  level,  the  greatest  bubble/bladder 
extensions  occur  when  cut-off  is  at  the  instant  of  minimum  bubble 
volume . * 


*  Footnote  next  page 


r 
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(3)  A  step  decrease  in  the  ambient  pressure  always  results 
in  an  extension  of  the  bubble/bladder.  If  such  a  drop  follows  a 
positive  step  where  cut-off  occurs  at  maximum  bubble  volume,  the 
extension  beyond  initial  at-rest  volume  is  due  solely  to  the  negative 
portion  of  the  pressure  signature;  if  not,  the  over-extension  is  due 
to  both  the  positive  step  and  the  underpressure  following  cut-off. 

The  effects  are  additive  in  the  non-linear  manner  described  in  this 
appendix . 


(Footnote  from  preceding  page) 

*For  a  simple  step  pressure  pulse  of  level  p  ,  the  energy  of  the  residual  bubble 
oscillation  is  given  by 


Y  =  Y_  -  Y.  =  (p  -  p.)  x  (V,  -  V  ) 

d  d  1  11  i 

where  V  is  the  bubble  volume  at  the  instant  of  cut-off.  Thus,  the  energy  of  the 


residual  oscillation  and,  consequently,  extension  of  the  bubble/bladder  beyond 
initial  at-rest  volume  take  on  greatest  values  for  V  =  V  and  tend  to  zero  as 


m 


V  approaches  V  . 
c  M 
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APPENDIX  B 

METHOD  FOR  CALCULATING  GAS  BLADDER  RESPONSE 
TO  EXPLOSION  PRESSURE  WAVE 

Figure  2.3.2  sketched  the  way  we  approximated  the  explosion 
pressure  signature  by  a  sequence  of  pressure  steps  in  order  to 
calculate  the  bladder  response  by  means  of  the  equations  developed 
in  Appendix  A.  Appendix  B  gives  the  details  of  the  procedure. 

Boundary  Condition  at  Pressure  Jumps.  Step  changes  in  the 
outside  water  pressure  from  pn  to  Pn+^  occurring  at  half-period 
intervals  change  the  equilibrium  pressure  PQ  of  the  oscillatinq 
bubble  flow  (equal  to  the  outside  pressure  p^) .  Since  the  oscillatory 
system  described  by  Equation  Al  has  finite  mass,  neither  A  nor  A  can 
change  impulsively,  i.e., 


A  A  =  0  (Bl) 

and 

A  A  =  0  (B2) 


for  all  pressure  jumps.  And,  thus  for  pressure  jumps  occurring  at 
extrema  where  A=0 ,  only  the  equilibrium  pressure  Pn  can  change. 

Let  A  be  the  maximum  or  minimum  bubble  radius  at  the 
cn 

time  of  the  pressure  jump  pn  to  Pn+^-  If  &c  is  greater  than  the 

n 

new  equilibrium  radius*  A  +^,  the  new  oscillation  begins  at  maximum 


*  _ 

The  equilibrium  radius  A  is  the  at-rest  radius  which  corresponds  to  the 
equilibrium  pressure  %  (=  p^,  the  outside  pressure).  A^  is  calculated 

from  p^  using  Equation  B6,  below. 


B-l 
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size;  and  if  Ac  is  smaller  than  An+^  the  new  oscillation  begins  at 
n 

minimum  size.  Or,  restated  we  have 


A  >  A  .  then  AMAX  .  =  A 
c  n+1  »  n+1  c 

n  n 


A  <  A  ,  then  AMIN 
c  n+1  >  n+1  =  A 

n  c 


Equilibrium  Bubble  Radius.  At  any  time  during  the  motion 
the  internal  gas  pressure  p  is  related  to  the  bubble  radius  by 


P(A)  =  p. 


Setting  P(A)  equal  to  p^  in  (B5)  yields  the  equilibrium  bubble 
radius , 


n  "Vn  / 


l/3y 


where  y  is  the  adiabatic  exponent  for  air  (=  1.40)  and  the  subscript 
"i"  refers  to  the  initial  ambient  "at  rest"  state. 

Equilibrium  Bubble  Period.  The  equilibrium  bubble  period 

T  (for  vanishingly  small  oscillations  about  the  equilibrium  radius) 

n  .  5 

was  given  by  Kennard. 


T  =  2ttA  /p/3Yp 
n  n  n 


where  p  is  the  water  density  =  1.940  slugs/ft“ . 
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Inserting  the  values  for  the  water  density  and  adiabatic  exponent, 
equation  B7  becomes 


T 

n 


milliseconds 


(B7a) 


where  An  is  in  inches  and  is  pounds  per  square  inch. 

Maximum  and  Minimum  Radius,  and  Period  of  Oscillation. 

With  increasing  amplitude  the  period  increases;  the  ratio  of  the 
periods,  T/T,  is  plotted  as  a  function  of  the  maximum  radius  AMAX 
in  Figure  Bl.  The  same  figure  also  shows  the  minimum  radius  AMIN 
plotted  as  a  function  AMAX.  Both  curves  were  calculated  by  the 
method  given  in  Appendix  A. 

Some  Details  of  the  Computation.  Consider  the  computation 
shown  in  Figure  2.3.2.  We  will  assume  that  the  pressure  signature 
begins  with  a  shock  front  followed  by  an  approximately  exponential 
decay.  To  get  the  average  outside  pressure  for  the  first  half-period 
of  the  motion  we  estimate  the  period  of  oscillation  by 


T  =  ex.  x  T* 


(B8) 


where  T*  is  the  equilibrium  period  calculated  using  the  pressure 
value  at  the  start  of  the  motion  (i.e.,  the  shockwave  peak  pressure  + 
initial  ambient)  and  is  an  arbitrary  constant.  (For  the  present 
computations  an  °<-value  of  1.25  was  used  for  the  first  half-cycle 
of  oscillation.) 

The  parameter  T*  is  calculated  from  Equation  B7  using  a  corresponding 
A*  calculated  from  Equation  B6 .  The  startinq  point  is  the  initial 
bubble  or  bladder  radius  A^  which  was  estimated  from  experimentally 
measured  values.  For  the  present  calculations  the  ratio  (A^)0/L 
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(Equations  3.1.7  and  3.1.8)  was  then  varied  to  obtain  the  best 
possible  correlation  between  the  calculated  and  experimentally 
observed  in juries--Section  3.2  and  3.3. 

Next,  we  read  or  calculate  the  outside  pressure  at  the  end 
of  the  first  half-cycle  of  motion  using  the  period  estimate  calculated 
by  (B8)  ;  and  then  calculate  p^ ,  the  constant  outside  pressure  for 
this  approximate  solution,  as  the  average  of  the  pressures  at  the 
beginning  and  end  of  the  half -cycle.  Using  pp  we  can  then  calculate 
the  parameters  for  the  first  half-cycle  of  motion;  and,  in  particular, 
obtain  a  considerably  improved  value  for  T.  Using  the  improved 
value  for  T  we  then  calculate  from  Equation  B8  a  better  value  for 
which  we  use  to  calculate  the  next  half-cycle  of  motion.* 

Surface  Cut-off.  Upon  arrival  of  the  reflected  wave  from 
the  water  surface  the  outside  water  pressure  suddenly  drops  to 
below  the  initial  ambient  level.  This  phenomenon  is  known  as 
"surface  cut-off".  When  cut-off  occurred  near  an  extremum — i.e., 
within  0.10  cycles  of  a  maximum  size,  or  within  0.08  cycles  of  a 
minimum--we  calculated  the  final  motion  at  the  lowered  outside 
pressure  as  starting  from  rest  from  the  nearby  extremum.  Interme¬ 
diate  cases  were  calculated  as  starting  from  rest  using  an  interme¬ 
diate  initial  bubble  size — based  on  the  parameter  (At/T)pc  which 
was  calculated  as  follows; 

If  A  +  AMAX  ,  then  (At/T)_,_  =  (TPOS  -  t  ) /T  (B9) 

cn-l  n  PC  cn-l  n 


*If  necessary,  we  can  use  this  improved  value  for  T  to  get  a  tetter  value  of  the 
average  pressure  p^.  So  far  this  has  not  been  necessary. 
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If  A  =  AMIN  ,  then  (At/T)  =  (TPOS  -  t  ) /T  +  1/2  (BIO) 

Cn-1  n  Cn-1  n 


The  parameter  (At/T)pc  locates  positive  phase  cut-off  in 
bubble  oscillation  cycles.  "Zero"  is  taken  as  the  last 


terms  of  bubble  oscillation  cycles.  "Zero"  is  taken  as  the  last 
expansion  occurring  during  the  positive  phase.  Using  (At/T)pCJthe 
starting  radius  for  the  final  oscillation  at  lowered  pressure  was 
then  selected  from  the  approximate  values  listed  in  Table  B-l. 


Phase  Shift  Correction.  The  effect  of  the  negative 

pressure  phase  is  often  limited  by  its  short  duration  relative  to 

the  bladder  period  of  oscillation.  Thus,  it  is  also  necessary  to 

account  for  the  phase  shift  which  occurs  upon  surface  cut-off.  This 

was  done  by  adding  a  time  correction  to  the  duration.  At  of  the 

neg, 

negative  pressure  phase.  Using  the  parameter  (At/T)pc  computed 
from  Equation  11  or  12,  this  correction  was  approximated  as  follows: 


If  0  £(At/T)pc  _<  0.1,  then  Phase  Corr  =  0  (Bll) 


If  0.1  <_  (At/T)pc  _<  0.5;then  Phase  Corr  =  -(t  -  TPOS)  (B12) 

n 


If  0.5  <_  (At/T)  <_  1.0;then  Phase  Corr  =  TPOS  -  t 


(B13) 


where  t  is  the  time  at  which  AMIN  would  have  occurred  had  not 
cn  n 

surface  cut-off  intervened.  (Since  this  subsequent  minimum  depends 
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TABLE  B-l 

APPROXIMATE  BOUNDARY  CONDITION  FOR  USE  WHEN  PRESSURE  JUMPS 

OCCUR  BETWEEN  EXTREMA 


START  NEXT 
HALF-CYCLE 
WITH: _ 

£0.10  AMAX 

<  0.25  (AMAX  +  A) /2 

£0.32  A 

<0.42  (A  +  AMINJ/2 
£0.58  AMIN 

<  0.68  (AMIN  +  A)/2 

£0.75  A 

<0.90  (A  +  AMAX) /2 

<  1.00  AMAX 
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primarily  on  the  kinetic  energy  at  cut-off  and  not  on  the  outside 
pressure,  this  tc~value  calculated  using  the  higher  pressure  should 
still  be  a  good  approximation  to  the  time  of  occurrence  of  this 
minimum. ) 


Final  Oscillation  at  Negative  Pressure.  In  calculating 
the  final  bubble  expansion  achieved  during  the  negative  pressure 
phase  we  only  consider  those  cases  in  which  the  final  oscillation 
starts  with  a  minimum  radius — only  they  are  of  interest.  In  order 
to  estimate  this  final  expansion  we  calculate  the  parameter 


(At/T)NP  =  (DTNEG  +  Phase  Corr)  /  Tne(?  +  1/2  (B14) 


where  DTNEG  is  the  duration  of  the  negative  phase — and  again  refer 
to  Table  B-l. 

The  parameter  (At/T)Np  locates  the  end  of  the  negative 
phase  in  terms  of  bubble  oscillation  cycles.  "Zero"  is  taken  as 

the  last  expansion  occurring  during  the  positive  phase. 

If  (At/T)Np  <_  0.58,  we  conclude  the  negative  phase  is  too  short  for 

any  significant  expansion  to  take  place.  Otherwise,  we  calculate 
the  maximum  size  as  indicated  in  the  table.  And,  if  (At/T)  p  •>  0.90, 

the  maximum  bubble  size  is  taken  equal  to  AMAX.  In  any  event,  if 
(At/T)  ^jp  <  0.90,  w«  also  calculate  the  final  oscillation  that  would 
occur  if  surface  cut-off  had  returned  the  outside  pressure  directly 
to  ambient.  If  this  results  in  a  greater  maximum  expansion,  we 
take  this  value  as  the  maximum  size  achieved  during  final  expansion. 

Damping .  Equation  Al  describes  undamped  radial  pulsation 
of  an  ideal  frictionless  fluid.  While  it  seems  unnecessary  to  make 
a  detailed  study  of  energy  dissipation  in  the  fishes'  oscillating 
swim  bladder,  it  does  seem  prudent  to  at  least  approximately  account 
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for  the  fact  that  such  dissipation  must  occur.  In  this  study  we 
did  so  by  withdrawing  a  fraction  B  of  the  energy  of  the  oscillatory 
motion  at  each  extremum. 

The  total  energy  at  an  extremum  can  be  calculated  from 
Equation  A1  by  dropping  the  first  term  and  substituting  from  (A2) 
and  (B5)  to  get 


v  -  ,  ^i  (  ^VMm^n  \ 

Yn  '  pn  v^T  (-V—  ) 


(VMm>n 


4  3 

where  is  the  bubble  volume,  y11  n 


Rewriting  (B5)  as 


p (V)  =  (v/vi) 


(B15) 


(B16) 


and  noting  that  pn  =  p(vn)  =  p^(Vn/V^)  ^  and  substituting 
into  (B15)  we  get 


Y  =  ^  p  v 

n  y-1  vn 


(B17) 


for  the  energy  of  the  non-oscillating  equilibrium  bubble. 


Subtracting  Yn  from  Yn  we  get  the  bubble  energy  of 


oscillation , 


(B18) 
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Finally  combining  (B15) ,  (B17)  and  (B18)  and  again  noting  that 

Pn  =  ( V/V^ )  Y  we  get  the  dimensionless  equation, 

V 

,Mm  _  (B19 ) 

V 

for  the  bubble  energy  of  oscillation. 

At  each  extremum  the  bubble  volume  (or  radius)  was  changed 
to  yield  a  new  energy  of  oscillation. 


'  =  Xzl  +  —  (  —  \ 


-y 


Y"  =  (1  -6)  Y' 
For  this  study  we  took 


(B20 ) 


6  =  0.30 


(B21) 


Considering  the  precision  of  our  present  experimental  data,  the 
precise  value  chosen  for  g  does  not  appear  to  be  critical. 
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APPENDIX  C 

APPROXIMATE  METHOD  FOR  CALCULATING  THE  PRESSURE-TIME  SIGNATURE 

Figure  3.1.1  shows  the  idealized  form  of  pressure-time 
signature  used  for  the  bladder  response  computations  of  this  report. 
For  planning  the  1975  series  of  explosion  tests  and  for  the  examples 
presented  in  Section  4  of  this  report  we  needed  to  predict  for 
arbitrary  explosion  geometries  the  variables,  PMAX,  0,  TPOS,  PNEG, 
DTNEG ,  shown  in  Figure  3.1.1.  This  Appendix  tells  how  this  was  done. 


PRESSURE  SIGNATURE  BEFORE  SURFACE  CUT-OFF.  The  initial 
positive  portion  of  the  signature  is  the  direct  wave  from  the  under¬ 
water  explosion.  The  sudden  negative  excursion  is  caused  by  the 
arrival  of  the  reflected  rarefaction  wave  from  the  water  surface. 
This  phenomenon  is  known  as  "surface  cut-off".  We  calculate  the 
time,  TPOS,  of  this  arrival  relative  to  the  direct  arrival  by 


TPOS 


(Cl) 


where  R  is  the  slant  range  from  the  charge,  R^  is  the  slant  range 
from  the  image  of  the  charge  mirrored  in  the  water  surface,  and  c0 
is  the  sound  speed  in  the  water.  See  Figure  Cl. 

We  calculate  the  free-field  portion  of  the  signature  using 
a  modification  of  the  empirical  Shockwave  Similitude  Equations.  The 
peak  pressure  PMAX  and  the  decay  constant  0  are  calculated  by  the 
usual  similitude  relations 
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PMAX 


(C2) 


(C3) 


where  W  is  the  explosive  weight  of  the  charge,  R  is  the  slant  range 
from  the  charge,  and  k,  a  ,  SL ,  and  B  are  empirical  constants  (see, 
e. g.  ,  Reference  7).  Values  for  the  constants,  k,  and  l  and  B  used 
for  the  calculations  of  this  report  can  be  obtained  from  Eguations 
3.1.1  and  3.1.3. 

The  instantaneous  free-field  (or  direct  arrival)  pressure, 
PD,  was  calculated  using  two  separate  exponential  seqments  joined 
at  t  =  1 . 8  0 


PD(t)  = 


PMAX  e_t//° 

0.25  PMAX  e-t'/4-3  0 


(t  <  1.8  0) 
(t  >  1.8  0) 


(C4a ) 
(C4b) 


which  were  fitted  to  the  pressure  signatures  recorded  on  the  1973 
explosion  test  series. 

CALCULATION  OF  PNEG,  NEGATIVE  PRESSURE  FOLLOWING  SURFACE 
CUT-OFF .  PNEG  was  determined  by  a  calculation  along  the  surface 
reflected  ray  from  the  charge  to  the  fish.  If  the  surface  reflected 
ray  went  thru  the  region  of  bulk  cavitation,  PNEG  was  computed  from 
the  depths  at  which  the  ray  intersected  the  top  and  bottom  of  the 
cavitation.  Otherwise,  PNEG  was  computed  from  the  linear  super¬ 
position  of  the  direct  and  surface-reflected  pressure  waves. 
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Calculation  of  PNEG  Above/Inside/Under  Region  of  Cavitation. 
Let  the  depths,  and  yc,  locate  the  top  and  bottom  of  the  region 
of  cavitation.  Let  p. __=p . +PNEG  be  the  absolute  pressure  under  the 

A£5o  1 

influence  of  cavitation,  i.e.,  following  arrival  of  the  surface 
reflected  wave.  For  fish  located  above/inside/under  the  region  of 
cavitation,  PABS  is  determined  as  follows: 


Above : 


PabS  Ganges  linearly  with 
the  surface  to  PVAP  at  the 


depth  from  PATM  at 
top  of  the  cavitation. 


Inside:  p.D<,  =  PVAP 

Ado 

Under:  PABS  =  PVAP  +  pg(y-yc) 

where  PVAP  and  p  are  the  vapor  pressure  and  density  of  water,  and 
g  is  the  acceleration  of  gravity.  Since 


PNEG  =  Pabs  -  Pi 


the  corresponding  equations  for  PNEG  are  (see  Figure  C2) 


PNEG  = 


/-[PATM  +  pgyt  -  PVAP ] 0<y<yt 


-[PATM  +  pgy  -  PVAP]  y^<Y^Y(. 


\  -[PATM  +  pgyc  -  PVAP] 


y>y. 


(C5) 


(C6a) 


(C6b) 


(C6c) 
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y 


FIG.  C-2  VARIATION  OF  PNEG  ABOVE/INSIDE/UNDER  REGION  OF  BULK  CAVITATION 
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Calculation  of  Depth,  y^,  of  Bottom  of  Cavitation.  Given 

Equations  C6a,  b,  c  the  problem  of  determininq  PNEG  reduces  to 
determination  of  the  depths  y  and  y^  of  the  top  and  bottom  of  the 
cavitation.  To  determine  y^  we  plotted,  as  a  function  of  the  charqe 
weight  W,  values  of  y  calculated  from  PNEG-values  measured  on  the 

c  i 

1973  Chesapeake  Bay  Test  Program  using  Equation  C6c.  The  equation.* 

yc  =  4.4  W0,3  (C7) 

where  W  is  the  explosive  weight  in  pounds  pentolite  and  y  is  in 
feet  was  used  to  approximate  these  measurements. 

Calculation  of  Depth,  y  ,  of  Top  of  Cavitation.  The  top 

(beginning)  of  the  cavitation  region  was  taken  as  the  point  on  the 
reflected  ray  where  the  reflected  shock  first  lowered  the  absolute 
pressure  to  the  vapor  pressure  of  the  water.  To  calculate  y  we  use 
the  linear  superposition  of  the  direct  and  surface-reflected  pressure 
waves , 


PSUM^  PD^  +  PR^ 

The  direct  wave,  pD(t),  is  given  by  equations  C 2,  3,  4.  The 
reflected  wave,  pR(t),  can  also  be  calcuated  usinq  Equations  C2,  3,  4 
provided  one  makes  the  substitutions  "-k"  for  "k"  and  "R^ "  (slant 
range  from  image)  for  "R".  To  determine  yfc  we  then  solve  for  the 
point  on  the  reflected  ray  path  where 


•Equation  C7  gives  a  crude  empirical  fit  to  the  PNEG-values  measured  on  the  1073 
test  program.  It  was  used  out  of  necessity  to  make  kill  probability  calculations 
used  to  design  the  1975  test  program  but  has  not  been  re-evauluated  using  those 
and  other  available  PNEG-measurements . 
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p.  +  p_  (TPOS)  +  p  (TPOS )  =  PVAP  (C9) 

1  U  K 

Calculation  of  PNEG  when  Reflected  Ray  is  Beyond  Region  of 

Cavitation .  If  the  total  pressure,  p.  +  p  (TPOS)  +  p  (TPOS), 

1  D  R 

upon  arrival  of  the  surface  refection  at  depth,  yc  =  4.4  ' 3 

(Equation  C7) ,  was  greater  than  or  equal  to  PVAP,  the  reflected  ray 
was  considered  to  lie  beyond  the  region  of  cavitation.  When  this 
happened  PNl’G  was  approximated  by 


PNEG  =  |(pD(TPOS)  +  pR(TPOS)] 


(CIO) 


calculated  at  the  location  of  the  fish. 

CALCULATION  OF  DTNEG ,  DURATION  OF  NEGATIVE  PRESSURE,  PNEG. 
DTNEG  was  also  determined  by  a  calculation  along  the  surface- 
reflected  ray  from  the  charge  to  the  fish.  If  the  surface-reflected 
ray  went  thru  the  region  of  cavitation  and  the  condition. 


R/W 


1/3 


<  80 


(Cll) 


where 


R  =  slant  range  from  charge  to  fish  in  feet 
W  =  explosive  weight  (pentolite)  in  pounds 

was  true,  then  DTNEG  was  calculated  from  the  time-of-f light  of  the 
water  layer  on  top  of  the  cavitated  region.  Otherwise,  DTNEG  was 
computed  from  the  linear  superposition  of  the  direct  and  surface- 
reflected  pressure  waves. 
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Calculation  of  DTNEG  from  Time-of-Flight  of  Water  Layer. 

The  duration,  DTNEG,  of  the  negative  phase  corresponds  to  the 
duration  of  bulk  cavitation  in  the  neighborhood  of  the  fish.  For 
this  study  the  duration  of  cavitation  was  calculated  at  the  point  on 
the  reflected  ray  at  depth  yc  (closure  depth)  given  by  Equation  C7 . 
Figure  C3  shows  the  geometry  of  the  problem.  To  get  a  first  approxi¬ 
mation  to  the  duration  of  bulk  cavitation  at  this  point  we  used 
Walker  and  Gordon's  result  for  time  of  flight  of  a  water  layer  of 

thickness,  y  ,  decelerating  (falling  back)  due  to  gravity  and  atmo- 

.  •  C  11 
spheric  pressure 


TFLIGHT 


2  •  PMAX  •  e 
_ _ _ c  c _ 

pg (y  +  k)  +  PATM 
c 


(C12 ) 


where  PMAXc  and  ©c  are  the  shockwave  peak  pressure  and  decay  constant 

at  the  closure  point  (x  ,y  )  , 

c  c 

and 


k  = 


9  R 
c  z 


DOB 


and 


c0  =  sound  speed  in  water 

Rz  =  slant  range  from  charge  to  point  of  reflection  at 
water  surface 
DOB  =  depth  of  burst. 


11/  Walker,  R.  R. ,  and  J.  D.  Gordon,  1966,  "A  Study  of  the  Bulk  Cavitation  Caused 
by  Underwater  Explosions",  David  Taylor  Model  Basin  Report  1896,  pl8. 

Equation  5*  Note,  however,  that  we  used  an  empirical  equation  ( CT )  to  calculate 
the  closure  depth  rather  than  Walker  and  Gordon's  result. 


IMAGE 
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To  determine  DTNEG,  we  plotted  the  ratio  of  TFLIGHT  (calculated  by 

C12)  to  DTNEG-values  measured  on  the  1973  Chesapeake  Bay  Program*  , 

as  a  function  of  R  /AMAX,  where  R  is  the  slant  ranqe  from  charae  to 

c  c 

cavitation  closure  point  and  AMAX  is  the  maximum  radius  of  the 
explosion  bubble.  This  plot  is  represented  by  the  empirical  equation^ 

9 

TFLIGHT 

DTNEG  = - (Cl  3) 

4.04  -  1.13  In (R  /AMAX) 
c 


where 

r  w  -.1/3 

AMAX  =  12.7  - - -  (C14 ) 

LdOB  +  (33.9/14. 7) -PATmJ 

and,  AMAX  and  DOB  are  in  feet,  PATM  is  in  pounds  per  square  inch, 
and  W  is  the  explosive  weight  in  pounds  pentolite. 

Thus,  in  summary,  when  the  reflected  ray  intersected  the 
region  of  cavitation  and  Equation  Cll  was  true,  Equations  C12,  13,  14 
were  used  to  calculate  DTNEG. 

Calculation  of  DTNEG  from  Superposition  of  Direct  and 
Surface-Refelcted  Waves.  If  the  reflected  ray  did  not  intersect  the 
region  of  cavitation  or  Equation  Cll  was  false,  DTNEG  was  calculated 
from  the  linear  superposition  of  the  direct  and  reflected  pressure 
waves  by  means  of  an  empirical  equation  adjusted  to  the  1973 
Chesapeake  Bay  Test  data.  This  empirical  equation  was 


•This  empirical  fit  was  used  to  make  the  kill  probability  calculations  for  design¬ 
ing  the  1975  test  program.  It  has  not  been  re-evaluated  using  the  1975  test 
program  and  other  available  DTNEG-measurements . 
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DTNEG  =  p'^M(PNEG/10)  -  TPOS 


(Cl  5 


where  PS^M  is  the  inverse  of  Equation  C8,  TPOS  is  the  positive 
duration,  and  PNEG  is  calculated  by  the  appropriate  equation,  C6  or 
CIO.  Thus,  the  end  of  the  negative  phase  was  taken  as  the  point 
where  the  negative  pressure  calculated  by  Equation  C8  returned  to 
10%  of  the  value,  PNEG. 


C-ll/12 
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APPENDIX  D 

LOCATION  OF  FISH  CAGES  BY  SOUNP  RANGING 

On  four  of  the  six  shots  of  the  1975  tests,  water  currents 
caused  significant  deviation  of  the  charge  and  the  fish  cages  from 
their  intended  positions.  The  charge  was  suspended  on  one  line  and 
the  cages — with  a  pressure  gage  attached  to  each--whore  suspended  on 
another.  In  this  appendix,  we  consider  the  problem  of  determining 
the  deviations  of  the  charge  and  the  string  of  fish  cages  from  their 
intended  locations. 

The  separation,  x0,  of  the  two  lines  at  the  surface  was 
assumed  known.  Figure  D1  shows  the  geometry  of  the  problem.  The 
length  of  the  charge  support  line  and  the  lengths  of  line  between 
gages  were  also  assumed  known.  We  further  assumed  that  the  charge 
and  gages  lie  in  the  the  same  vertical  plane. 

We  solved  this  problem  using  the  measured  arrival  times 
of  the  direct  shockwave  and  its  reflected  tension  wave  from  the 
water  surface  at  the  respective  gage  locations.  Step  1  was  to  use 
the  relative  arrivals  of  the  direct  shock  to  determine  an  orientation 
for  the  gage  line  giving  a  set  of  constant  differences  between  the 
calculated  arrivals  and  the  measured  ones.  Step  2  was  to  use  these 
gage  positions  to  calculate  the  surface  cut-off  time — the  time 
interval  between  arrival  of  the  direct  shock  and  its  reflection — 
and  to  compare  these  values  with  the  corresponding  measured  values. 

If  there  existed  a  systematic  discrepancy  in  surface  cut-off  times 
not  accountable  to  variations  in  the  sound  speed  over  the  two  ray 
paths,  we  then  adjusted  the  angle,  0  ,  of  the  charge  support  line 
and  redid  steps  1  and  2.  The  calculations  were  redone  until  there 
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was  no  systematic  discrepancy  in  surface  sut-off  time  not  accountable 
to  variation  in  sound  speed  over  the  two  ray  paths. 

The  computations  of  steps  1  and  2  assumed  that  the  sound 
speed  did  not  vary  with  depth.  R.  S.  Price  of  this  center  has 
developed  a  high-speed  computer  program  which  takes  into  account  the 
sound  speed  variation  with  depth--but  assumes  that  the  gages  are  in 
a  straight  line.  For  the  computations  done  here,  these  two  effects — 
curvature  of  the  gage  line  and  sound  speed  variation — are  the  same 
order  of  magnitude.  In  these  computations,  errors  of  about  0.4 
milliseconds  are  caused  by  variations  in  sound  speed.  These 
correspond  to  gage  position  errors  of  about  2  feet.  For  our  purposes 
errors  of  this  magnitude  were  acceptable,  however,  significantly 
improved  precision  could  be  obtained  by  accounting  for  both  gage  line 
curvature  and  sound  speed  variation  in  the  calculations. 

As  a  final  check  on  the  computations  we  compared  the 

absolute  time  of  arrival  of  the  direct  shockwave  measured  on  the 

pressure  gage  records  to  that  calculated  for  the  shockwave  travelling 

between  our  computed  charge  and  gage  locations.  To  do  this  we 

calculated  the  time,  At  ,  between  the  electrical  firing  pulse 

meas  r 

and  direct  shockwave  arrival  as  follows 


At  _  =  At  ,  +  VEL.  CORR.  +  FIRING  DELAY 
meas  sound 


(Dl) 


where 


At 


sound 


transit  time  for  a  sound  wave 


VEL.  CORR.  =  correction  to  At  ,  to  account  for  detonation 

sound 

wave  velocity  inside  the  charge  and  shockwave 
velocity  in  the  water 

FIRING  DELAY  =  dwell  time  between  firing  pulse  and  initiation 
of  charge 
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The  firing  delay  depends  on  the  firing  circuit,  length  and 
type  of  firing  line,  and  the  detonator.  For  the  1975  tests  we 
estimated  from  a  comparable  test  setup  that 


FIRING  DELAY  =  +0.31  +  0.05  milliseconds 


(D2) 


The  correction,  VEL.  CORR.  ,  to  must  be  calculated. 

It  is  given  by  the  following  equation 


VEL.  CORR. 


+ 


where 


(D3) 


R0  =  radius  of  the  charge 

c0  =  sound  velocity  in  the  water 

U  =  shock  velocity  in  the  water 

R  =  radial  distance  to  gage 

D  =  detonation  velocity  in  the  charge 

The  first  term  inside  the  brackets  represents  the  contribution  due 
to  the  water  shock  while  the  second,  that  due  to  the  detonation  wave. 
Equation  D3  was  derived  from  the  sketch  shown  in  Figure  D-2 . 

The  integrand  in  Equation  D3  was  obtained  by  fitting 

results  from  a  hydrodynamic  code  calculation  for  pentolite  by 

12 

Sternberg  and  Walker  with  the  following  equations 


12/ 

—  Sternberg,  H.  M. ,  and  W.  A.  Walker,  1971,  "Calculated  Flow  and  Energy 
Distribution  Following  Underwater  Detonation  of  a  Pentolite  Sphere", 
Physics  of  Fluids,  September  1971 


FIG.  D-2  DETONATION  AND  SHOCKWAVE  VELOCITY  CORRECTIONS  TO  ARRIVAL  TIME 
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EXP  |8-84-]60-’  +  '  } 


EXP  | -3 . 264  -  1.098 


InB^j 


for  R/R0  <  25 


for  R,/R0  >  25* 


(D4a) 


(D4b) 


To  calculate  VEL.  CORR.  Equation  D3  must  be  integrated  all  the  way 
out  to  the  gage.  (This  correction  does  not  assume  a  constant  value 
beyond  some  arbitrary  point,  i.e.,  for  this  application  the  shock- 
wave  cannot  be  considered  to  be  like  a  sound  wave,  anywhere.) 

Table  D-l  summarizes  the  results  of  the  sound  ranging 
computations.  The  last  two  columns  give  the  comparisons  for  the 
time  of  arrival  for  the  direct  shockwave.  The  comparison  is  easily 
understood  by  subtracting  At  un<j  from  both  sides  of  Equation  D1 . 

The  maximum  discrepancy  is  0.4  milliseconds  (Shot  786)  or  a  spatial 
discrepancy  of  about  2  feet. 

Tables  D-2  thru  D-7  list  cage  and  gage  locations  for  the 
1975  test  series  which  where  calculated  by  the  method  of  this  appendix. 
The  values  listed  in  Tables  D-2  thru  D-7  were  used  for  the  present 
analysis  of  the  1973  and  1975  test  series  data.  The  distance  along 
the  support  wire  (Column  1)  is  the  nominal  or  desired  depth  of  the 
fish  cage  or  pressure  gage.  The  wire  angle  (Column  2)  is  the 
calculated  deviation  from  the  vertical  of  the  wire  segment  attached 
to  the  preceding  depth  coordinate.  The  cage  and  gage  coordinates 
(Columns  3  thru  6)  are  calculated  relative  to  the  point  on  the  water 
surface  directly  above  the  charge.  Sound  speeds  were  calculated 
from  salinity  and  water  temperature  at  charge  depth  using  Del  Grosso's 
equation. 1 ■  (All  the  computations  described  in  this  appendix  were 
done  using  an  HP-65  Programmable  Pocket  Calculator.) 

"This  straight  line  log-log  extrapolation  beyond  the  range  of  Sternberg  and  Walker's 
calculation  should  be  adequate  out  to  at  least  R/R0  =  1000. 

13. Del  Grosso,  V.  A.,  "New  Equation  for  the  Speed  of  Sound  In  Natural  Waters  (with 
Comparisons  to  Other  Equations)",  J.  Aeoust.  Soc.  Am.,  October  107)* 
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TABLE  D-2 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  782 

Charge:  70.4  lbs  pentolite 
Charge  wire  length:  30  feet 
Charge  wire  angle:  0  degrees 
Charge  depth:  30  feet 

Charge-to-gage  line  distance:  300  feet 
Sound  speed:  4850  ft/sec 
Water  depth*:  156  feet 


WIRE 

DISTANCE 

WIRE 

ANGLE 

CAGE  AND  GAGE  COORDINATE 

:s 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 

300 

0.0 

91.4 

0.0 

5 

0 

300 

5.0 

91.4 

1.5 

10 

0 

300 

10.0 

91.4 

3.0 

15 

0 

300 

15.0 

91.4 

4.6 

40 

0 

300 

40.0 

91.4 

12.2 

45 

0 

300 

45.0 

91.4 

13.7 

50 

0 

300 

50.0 

91.4 

15.2 

55 

0 

300 

55.0 

91.4 

16.8 

57.5 

0 

300 

57.5 

91.4 

17.5 

77.5 

0 

300 

77.5 

91.4 

2  3.6 

87.5 

0 

300 

87.5 

91.4 

26.7 

97.5 

0 

300 

97.5 

91.4 

29.7 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D  -  3 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  783 

Charge:  70.2  lbs  pentolite 

Charge  wire  length:  30  feet 
Charge  wire  angle:  0  degrees 

Charge  depth:  30  feet 

Charge- to-gage  line  distance:  200  feet 
Sound  speed:  4857  ft/sec 
Water  depth*:  152  feet 


WIRE 

DISTANCE 

WIRE 

ANGLE 

CAGE  AND  GAGE  COORDINATES 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

( DEG . ) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 

— 

200.0 

0.0 

61.0 

0.0 

5 

-10.5 

199.1 

4.9 

60.7 

1.5 

10 

-10.5 

198.2 

9.3 

60.4 

3.0 

15 

-10.5 

197 . 3 

14.7 

60.1 

4.5 

40 

-10.5 

192.7 

39.3 

58.7 

12.0 

45 

-10.5 

191.8 

44.2 

58.5 

13.5 

50 

-10.5 

190.9 

49.2 

58.2 

15.0 

55 

-10.5 

190.0 

54.1 

57.9 

16.5 

60 

-10.5 

189.1 

59.0 

57.6 

18.0 

80 

-10 . 5 

185.4 

78.7 

56.5 

24 . 0 

90 

-  8 

184.0 

88.6 

56.1 

27.0 

100 

-  5 

183.2 

— 

98.5 

55.8 

30.0 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D-4 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  784 

Charge:  71.6  lbs  pentolite 

Charge  wire  length:  30  feet 
Charge  wire  angle:  0  degrees 

Charge  depth:  30  feet 


Charge-to-gage  line  distance:  300  feet 
Sound  speed:  4866  ft/sec 
Water  depth*:  135  feet 


WIRE 

DISTANCE 

WIRE 

ANGLE 

CAGE  AND  GAGE  COORDINATES 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

(DEG. ) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 

— 

300.0 

0.0 

91.4 

0.0 

5 

11 

301.0 

4.9 

91.7 

1.5 

10 

11 

301.9 

9.8 

92.0 

3.0 

20 

10 

303.6 

19.7 

92.5 

6.0 

30 

10 

305.4 

29.5 

93.1 

9.0 

40 

9 

306.9 

39.4 

93.5 

12.0 

45 

8 

307.6 

44.3 

93.8 

13.5 

50 

6 

308.2 

49.3 

93.9 

15.0 

55 

4 

308.5 

54 . 3 

94.0 

16.6 

60 

3 

308.8 

59.3 

94.1 

18.1 

80 

2 

309.5 

79.3 

94.3 

24.2 

90 

1 

309.6 

89.3 

94.4 

27.2 

100 

1 

309.8 

99.3 

94.4 

30.3 

Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D-5 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  785 

Charge:  1.25  lbs  pentolite 
Charge  wire  length:  30  feet 
Charge  wire  anqle:  -5  degrees 
Charge  depth: 29.9  feet 

Charge- to-gage  line  distance:  40  feet 
Sound  speed:  4864  ft/sec 
Water  depth*:  160  feet 


WIRE 

DISTANCE 

WIRE 

ANGLE 

CAGE  AND  GAGE  COORDINATES 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 

42.6 

0.0 

13.0 

0.0 

5 

41.7 

4.9 

12 . 7 

1.5 

10 

40.7 

9.8 

12.4 

3.0 

20 

-ii 

39.3 

19.7 

12.0 

6.0 

30 

-  7 

38.1 

29.6 

11.6 

9.0 

40 

-  7 

37.1 

39.6 

11.3 

12.1 

50 

-  5 

36.0 

49.5 

11.0 

15.1 

60 

-  4 

35.0 

59.5 

10.7 

18.1 

70 

-  2 

33.9 

69.4 

10.3 

21.2 

80 

+  2 

34.1 

79.4 

10.4 

24.2 

90 

+  5 

34.4 

89.4 

10.5 

27.2 

100 

+  5 

35.1 

99.4 

10.7 

30.3 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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TABLE  D-6 


CORRECTED  CAG 


AND  GAGE  LOCATIONS 


SHOT  786 


Charge:  1.26  lbs  pentolite 


Charge  wire  angle:  -10  degrees 


Charge  depth:  29. 5  feet 


Sound  speed:  4886  ft/sec 


Water  depth 


WIRE 

ANGLE 


WIRE 

DISTANCE 


FEET) 


♦Calculated  from  times  of  arrival  of  bottom  reflection 


CAGE  ANC 

| 

)  GAGE  COORDINATES 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

“7  FEET)' 

(METERS) 

(METERS) 

45.2 

0.0 

13.8 

0.0 

42 . 7 

4 . 3 

13.0 

1.3 

40.2 

8.7 

12.3 

2.7 

36.1 

it.  8 

11.0 

5.4 

32.9 

27.3 

10.0 

8.3 

30.3 

36.9 

9.2 

11.2 

28.2 

46 . 7 

8.6 

14.2 

27.4 

51.6 

8.4 

15.7 

26.5 

56 . 5 

8.1 

17.2 

25.8 

61.5 

7.9 

18.7 

25.3 

66.5 

7.7 

20.3 

1  24.9 

71.5 

7.6 

21.8 

24. T 

76.4 

7.5 

23.3 
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TABLE  D  -  7 

CORRECTED  CAGE  AND  GAGE  LOCATIONS 
SHOT  787 


Charge:  72.2  lbs  pentolite 

Charge  wire  length:  Id  feet 
Charge  wire  angle:  0  degrees 
Charge  depth:  10  feet 

Charge-to-gage  line  distance:  300  feet 
Sound  speed:  4922  ft/sec 
Water  depth*:  148  feet 


— 

WIRE 

DISTANCE 

— 

WIRE 

ANGLE 

CAGE  AND  GAGE  COORDINATES 

HOR. 

RANGE 

DEPTH 

HOR. 

RANGE 

DEPTH 

(FEET) 

■mai 

(FEET) 

(FEET) 

(METERS) 

(METERS) 

0 

300 

0 

91.4 

0.0 

5 

0 

300 

5 

91.4 

1.5 

10 

0 

300 

10 

91.4 

3.0 

40 

0 

300 

40 

91.4 

12.2 

50 

0 

300 

50 

91.4 

15.2 

55 

0 

300 

55 

91.4 

16.8 

60 

0 

300 

60 

91.4 

18.3 

70 

0 

300 

70 

91.4 

21.3 

80 

0 

300 

80 

91.4 

24.4 

90 

0 

300 

90 

91.4 

27.4 

100 

0 

300 

100 

91.4 

30.5 

♦Calculated  from  times  of  arrival  of  bottom  reflection. 
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APPENDIX  E 

NOTE  ON  FISH  CLOSE  TO  THE  WATER  SURFACE 


According  to  the  bladder  oscillation  model  for  explosion 
injury,  as  the  fishes'  depth  tends  to  zero  near  the  water  surface 
the  kill  probability  should  also  tend  to  zero  reqardless  of  how  high 
the  incident  pressure.  Eventually,  however,  for  fish  at  the  water 
surface  and  sufficiently  close  to  the  charge,  the  bladder  oscillation 
model  must  break  down  and  a  different  mechanism  for  injury  must  take 
over.  Apparently,  this  happens  so  close  to  the  charge,  and  to  the 
water  surface,  to  be  of  little  practical  importance,  since  results 
reported  by  the  Lovelace  Foundation  ,4  show  no  evidence  of  such  a 
transition  for  149  gm  Carp  (length  =  21  cm)  at  a  depth  of  5  cm 
subjected  to  5.58  megapascals  (810  psi)  incident  pressure. 

A  possible  close-in  non-bladder  injurv  mechanism  is  tissue 
damage  due  directly  to  compression  by  the  incident  shockwave.  If 
this  is  the  case,  tissue  damage  probablv  also  occurs  due  to  the 
pressure  wave  (generally  an  order  of  magnitude  areater)  emitted  by 
the  oscillating  swim  bladder.  Assuming  that  this  is  the  meaning  of 
the  ratio,  AMIN/A^ ,  in  the  damage  parameter  Z,  we  rewrite  Eauation 
3.1.3  in  terms  of  the  adiabatic  compression  ratio,  PMIN/p0 ,  as 


X  =  -100  In 


(El) 


where  PMIN  is  the  maximum  value  of  the  oscillating  bladder  pressure 
(corresponding  to  radius  AMIN  at  the  first  compression) ,  p^  is  the 
initial  ambient  pressure,  and  y, the  adiabatic  exponent  for  air. 
Using  El  we  then  rewrite  3.1.5  as 
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Z  =  i--  In  PMIN/p .  +  100  In  AMAX/A.  (E2) 

JY  l  l 

Note  that  Equation  E2  is  equivalent  to  3.1.5.  It  is  just  our 
point  of  view  that  has  changed. 

Now  let  us  suppose,  for  fish  very  close’to  the  water 
surface,  that  the  shockwave  peak  overpressure  PMAX  is  the  only 
tissue  damaging  pressure.  The  damage  parameter  by  this  mechanism 
equivalent  to  Equation  E2  is  then  given  by 


''PMAX 


( E  3 ) 


where  PMIN  is  replaced  by  PMAX  +  p .  and  the  second  term  drops  out. 
Note  that  must  always  be  smaller  than  Z  calculated  by 

Equation  3.1.5  except  for  those  cases — the  ones  of  interest  here-- 
where  the  pressure  wave  is  of  such  short  duration  that  the  transient 
response  of  the  swimbladder  is  suppressed*.  Thus,  when  in  doubt 
which  damage  parameter  to  use--Zp^AX  (Equation  E3)  or 
Z  (Equation  3 . 1 . 5) --calculate  both  parameters,  and  use  the  greater 
va 1 ue . 

Finally,  we  estimate  the  overpressure  level,  PMAX, 
associated  with  the  assumption  that  all  the  damage  corresponding 
to  the  damage  parameter  is  due  directly  to  tissue  compression 

by  the  incident  shockwave.  Solving  E3  for  PMAX  we  get 
* - - - 

Since  zpj^x  is  equivalent  to  the  damage  parameter  X 

(Equation  3.1.3)  calculated  from  the  non-osci 1 1 atory  steady-state 
response  of  an  infinitely  small  bubble  to  an  overpressure  PMAX 
of  finite  rise  time. 
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PMAX  =  p.  EXP  [  (3YZpMAX/100)  -  l] 


Setting  p^  =  1  atmosphere  (1.014  x  10J  pascals)  and 
ZPMAX  =  (50%  kill  value--Equation  3.2.1)  we  get 

PMAX  =  19.2  MPa  (=2790  psi) 

(approximately  24  charge  diameters) 


for  the  overpressure  level  PMAX  corresponding  to  50%  kill 
by  this  mechanism  to  fish  which  are  very  close  to  the  water 
surface . 


E-3/E-4 


i 


k 


-n 
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